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ABSTRACT 


Thi  s 

investigation 

is  part 

of 

the 

current  research 

effort  in 

geotechni ca 1 

engineer i ng 

to 

integrate 

the 

deformat  ion 

and  stabi 1 i ty 

ana  1 yses 

of 

soi  1 

structures . 

Wi  th 

the  aid  of  recent  advances  made  in  the  mathematical 
modelling  of  soil  behavior,  the  engineering  profession  is 
now  a  step  closer  to  making  successful  predictions  for  the 
movements  of  soil  masses  at  all  stages  of  loading  up  to 
failure.  The  study  presented  here  is  an  attempt  to  evaluate 
the  performance  of  Lade's  work  hardening  stress-strai n 
relation  in  an  integrated  analysis.  Previous  applications  of 
the  model  carried  out  by  others  have  not  been  fully 
satisfactory  despite  several  corrections  and  proposed 
modi fi cat  ions . 

To  begin  with,  the  original  formulation  of  the  model 
and  subsequent  modifications  are  re-evaluated  as  part  of  the 
present  study.  It  is  found  that  the  difficulties  encountered 
in  the  past  are  due  to  the  nature  of  the  hardening  law  of 
the  model.  For  a  general  three  dimensional  problem,  a  unique 
relation  between  stress  level  and  plastic  work  does  not 
exist.  It  is  shown  here  that,  by  making  a  suitable 
assumption,  measured  stress-strain  behavior  of  cohesionless 
soils  can  be  successfully  predicted.  At  the  same  time,  the 
problems  that  are  created  in  finite  element  analysis  can  be 
avoided . 

For  further  evaluation  of  the  model,  a  finite  element 
program  has  been  developed  to  predict  the  results  of  passive 
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earth  pressure  tests  carried  out  by  Wong.  The 
load-deformation  response,  as  well  as  the  distribution  of 
normal  stresses  acting  on  a  fully  instrumented  wall  has  been 
predicted  with  sufficient  accuracy  for  all  practical 
purposes.  These  results  have  shown  that  Lade's  work 
hardening  model  works  very  well  for  the  analyses  of  complex 
boundary  value  problems  at  all  stages  of  loading.  When 
compared  with  other  models,  such  as  the  hyperbolic 
stress-strain  law  and  the  stress-di 1 atancy  model,  Lade's 
model  gave  much  better  results  for  the  same  wall  behavior. 

It  is  believed  that  the  major  factor  in  successful 
predictions  of  Lade's  work  hardening  stress-strain  relation 
is  the  reasonably  accurate  modelling  of  basic 
characteristics  of  soil  behavior. 
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CHAPTER  1 


INTRODUCTION 


1 . 1  GENERAL 

In  classical  soil  mechanics,  the  problems  of 
deformation  and  stability  are  treated  separately.  Combined 
with  experience,  these  methods  provide  satisfactory  results 
for  most  geotechnical  engineering  problems.  However,  the 
response  of  a  soil  mass  to  loading  is  one  continuous 
process,  and  the  classical  approach  is  therefore  artificial. 
In  some  deformation  problems,  it  is  quite  possible  to  have 
local  yielding  even  at  high  factors  of  safety.  Similarly, 
when  collapse  takes  place,  not  all  elements  of  a  soil  mass 
are  in  a  state  of  failure;  some  regions  may  still  be 
undergoing  loading  or  even  unloading.  Considering  that  the 
aim  of  a  stability  analysis,  at  least  for  some  problems,  is 
to  find  the  most  economical  way  of  avoiding  failure,  the 
collapse  load  needs  to  be  predicted  as  accurately  as 
possible.  Further,  the  mechanism  of  failure  is  mainly  a 
function  of  the  stress  distribution  which  can  only  be  found, 
in  general,  by  using  a  deformation  analysis  capable  of 
simulating  all  loading  stages  prior  to  instability. 

One  of  the  principal  requirements  for  an  integrated 
analysis  of  deformations  and  failure  is  a  constitutive 
relationship  suitable  for  modelling  the  st ress-strai n 
behavior  of  soil  correctly  up  to  and  beyond  the  peak 
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strength.  This  challenge  has  been  taken  by  a  growing  number 
of  researchers  since  the  first  establishment  of  the  Cam-clay 
model  ( Roscoe  et  al,  1958)  which  marked  the  beginning  of  a 
modern  approach  to  modelling  and  fundamental  understanding 
of  the  mechanical  behavior  of  soils. 

Soil  behavior  under  controlled  laboratory  conditions 
has  been  investigated  in  great  detail  and  is  well  documented 
in  the  literature.  Even  under  these  ideal  conditions  soil 
behavior  is  extremely  complicated  in  comparison  to  other 
engineering  materials.  Consequently,  constitutive  laws  which 
aim  to  model  most  aspects  of  soi 1  behavior  tend  to  become 
highly  complex.  Considering  that  the  behavior  of  natural 
soils  is  affected  by  many  other  factors  such  as  anisotropy, 
non -homogene i ty ,  time  and  temperature  effects,  the  task  of 
developing  mathematical  models  becomes  very  demanding  if 
every  aspect  of  soil  behavior  is  to  be  accounted  for. 

While  the  importance  of  developing  constitutive  models 
is  indicated  and  the  complexity  of  the  developer's  work  is 
acknowledged,  the  size  of  the  remaining  work  required  to 
complete  the  analysis  of  an  engineering  structure  should  not 
be  underestimated.  The  analysis  of  any  continuum  mechanics 
problem  requires  that  the  field  equations  must  be  solved 
subject  to  the  appropriate  boundary  conditions.  Also,  the 
effects  of  the  initial  state  of  stress  and  the  stress 
history  of  in-si tu  soils  have  to  be  considered  in  a 
realistic  analysis.  Finite  element  techniques  employing 
elastic  or  nonlinear  elastic  models  have  been  used  widely  in 
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the  past  for  the  analysis  of  important  soil  structures,  and 
considerable  experience  has  already  been  accumulated  in  this 
area.  But  much  more  effort  is  required  to  use  the 
complicated  constitutive  models  in  finite  element  analysis 
with  the  same  level  of  confidence  gained  so  far  for  linear 
and  nonlinear  elastic  solutions.  It  should  be  noted  that 
nonlinear  analysis  by  the  finite  element  method  is  still  an 
ongoing  research  area  requiring  a  high  level  of  expertise. 

In  the  following,  the  stress-strain  laws  which  attempt 
to  model  the  essential  features  of  soil  behavior  such  as 
nonlinearity,  inelasticity,  shear  dilatancy  and  path 
dependency  are  referred  to  as  "non-classical"  models.  This 
term  is  used  in  this  thesis  not  only  to  point  out  the 
remarkable  abilities  of  these  models  to  predict  the 
stress-strain  behavior  of  soils,  but  also  to  suggest  that 
their  application  to  the  analysis  of  soil  structures  may  not 
be  straight  forward  but  full  of  difficulties. 

It  would  seem  that  sufficiently  accurate  answers  can  be 
found  to  any  geotechnical  problem  once  the  finite  element 
programs  for  non-classical  stress-strain  laws  are  developed. 
However,  one  should  be  cautious  about  reaching  this 
conclusion  too  quickly.  Not  too  long  ago  Wroth  (1976)  stated 
that: 


"It  might  be  thougt  that 
in  foundation  engineering 


any  boundary  value  problem 
such  as  the  behavior  of 
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the  ground  around  an  excavation  could  now  be 
completely  and  accurately  solved  by  a  deformation 
analysis  using  the  finite  element  method.  But  it  is 
my  contention  that  is  not  so;  I  intend  to  show  that 
the  finite  element  method  as  developed  at  present  is 
unable  to  represent  in  an  adequate  manner  the 
discontinuous  nature  of  the  behavior  of  soils  at 
f a i 1 ure ..." 


In  addition  to  our  inability  to  introduce  the 
development  of  failure  surfaces  into  finite  element 
analysis,  geotechnical  engineering  practice  is  also  troubled 
in  determining  accurately  the  state  of  stress  and  the  stress 
history  of  in-situ  soils.  This  vital  information  is  the 
starting  point  to  any  analysis. 

Knowing  that  no  single  constitutive  law  describes 
adequately  all  features  of  soil  behavior,  the  analyst  has  to 
make  a  choice  of  the  most  appropriate  model  which  inevitably 
becomes  a  compromise.  Therefore  the  finite  element  method  of 
analysis  will  be  useful  to  engineering  practice  only  if  the 
limitations  of  particular  soil  model  and  solution  method  are 
recognized . 

It  should  also  be  noted  that  the  appreciation  of  all 
the  refinements  of  recently  developed  constitutive  laws  and 
the  understanding  of  the  limitations  of  finite  element 
analysis  are  not  the  only  issues  the  user  or  practicing 
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engineer  has  to  deal  with.  As  pointed  out  by  Morgenstern 
( 1975) , 


"Stress-strain  relations  for  soils,  either 
theoretical  or  empirical,  constitute  only  one 
element  in  the  total  design  procedure  that  involves 
at  the  very  least  calculations,  judgement,  quality 
control  and  monitoring  of  per formance . " 


In  relation  to  the  application  of  stress-strai n  laws 
there  are  several  important  points  to  be  considered.  Some  of 
them  are  discussed  below. 

1.  So  far,  there  has  not  been  any  comparative  study 
available  as  a  guideline  to  the  users  to  show  the 
capabilities  and  limitations  of  non-classical  models. 
Not  only  predictions  of  laboratory  test  results  but  also 
analyses  .of  the  same  engineering  structures  by  all 
models  considered  are  necessary  before  the  true  value  of 
each  model  can  be  assessed.  Therefore  it  is  difficult  to 
decide  among  a  dozen  or  more  different  constitutive 
models  which  one  to  choose  for  a  particular  problem  in 
hand,  unless  information  on  comparative  performance, 
cost,  availability  etc  has  been  provided. 

2.  For  each  model,  testing  techniques  required  and 

difficulties  involved  to  obtain  soil  parameters  should 
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be  the  first  questions  to  be  asked.  While  some  models 
require  several  different  types  of  tests,  for  others 
conventional  triaxial  tests  are  sufficient  for  deriving 
soil  parameters. 

3.  Due  to  the  difficulties  associated  with  representative 

sampling  and  sample  disturbance,  laboratory  tests  often 
are  of  limited  value  (Lambe,  1973).  Therefore  the 
existence  of  a  correlation  between  the  results  of 
present  day  in-si tu  testing  techniques  and  the 

parameters  of  a  stress-strai n  model  has  practical 
importance . 

4.  Some  of  the  non-classical  constitutive  models  depend 

heavily  on  curve  fitting  techniques  to  obtain  soil 
parameters  from  laboratory  test  results.  However,  due  to 
the  nature  of  soils,  it  is  difficult  to  duplicate  the 
test  results  for  the  same  type  of  soil  even  under  ideal 
laboratory  conditions.  Therefore  the  sensitivity  of  each 
soil  parameter  to  any  slight  variation  in  soil 

properties  and  also  to  the  accuracy  of  the  soil  testing 
technique  should  be  fully  investigated.  For  those  that 
vary  significantly  and  have  a  considerable  effect  on  the 
results  of  an  analysis  have  to  be  determined  and  the 
analysis  has  to  be  carried  out  accordingly. 

5.  It  is  not  always  necessary  to  use  the  most  complete 
stress-strain  model  in  every  analytical  problem.  The 
significance  of  modelling  certain  soil  characteristics 
depends  on  the  type  of  boundary  value  problem  and  the 
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nature  of  the  in-si tu  soils.  For  example  in  some 
problems  the  time  factor  becomes  important.  For  some 
others  cyclic  loading  does  not  exist.  A  constitutive  law 
capable  of  modelling  strain  softening  may  be  attractive 
but  a  work  hardening  model  may  well  be  satisfactory  for 
some  problems.  The  use  of  a  non-associ ated  flow  rule 
does  not  necessarily  produce  better  results  over  an 
associated  flow  rule  for  all  problems.  Similarly,  a 
failure  criterion  which  has  a  curved  envelope  in  the 
triaxial  plane  and  can  also  account  for  the  effect  of 
intermediate  principal  stress  may  not  be  significantly 
superior  to  the  Mohr-Coulomb  failure  criterion.  So  far 
there  has  not  been  any  systematic  study  of  real  soil 
structures  demonstrating  the  advantages  of  the  effort 
involved  to  include  the  refinements  mentioned  above  in 
an  analysis.  A  parametric  study  using  the  type  and  the 
size  of  soil  structure,  the  boundary  conditions,  and  the 
kind  and  the  state  of  soil  as  variables  would  be 
valuable  for  engineering  practice. 

6.  Although  some  models  are  claimed  to  account  for  most 
characteristics  of  soil  behavior,  they  are  developed 
mostly  on  the  basis  of  laboratory  tests  on  kaolinite  or 
uniform  sand  (“university  materials").  The  application 
of  such  models  to  the  analysis  of  soil  structures  where 
natural  materials  exist  may  be  a  good  test  for  those 
mode  1 s . 

7.  Once  the  type  of  model  to  be  used  has  been  decided,  the 
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availability  of  the  finite  element  program  specifically 
written  for  the  model  has  to  be  ascertained.  At  this 
point,  the  difficulties  involved  in  its  use  have  to  be 
determined.  If  modifications  are  required  or  problems 
are  encountered  during  analysis,  it  is  usually  necessary 
to  consult  the  developer  of  the  program  or  even  the 
developer  of  the  model  itself.  These  processes  may 
greatly  complicate  the  application  of  a  model. 

8.  The  cost  and  the  time  required  to  run  a  finite  element 
program  depend  on  the  complexity  of  the  model  adopted 
and  the  refinement  necessary  to  solve  that  particular 
problem.  This  can  be  a  decisive  factor  in  choosing  a 
soil  model  for  analysis  depending  on  the  size  and  the 
importance  of  the  project. 

9.  Sometimes  the  accuracy  of  the  results  of  a  finite 
element  program  becomes  questionable.  Depending  on  the 
numerical  technique  used  to  solve  the  nonlinear  problem, 
there  are  always  approximations  introduced  into  the 
ana  lysis. 

10.  The  reliability  of  the  finite  element  program  at  hand 
should  be  checked  carefully.  For  any  situation  different 
from  a  previous  successful  use,  the  program  should  be 
treated  as  suspect  until  proven  otherwise.  Every  program 
has  to  go  through  a  testing  period  before  its  results 
are  relied  upon . 

The  points  made  in  this  introductory  section  form  the 

general  framework  of  present  study.  An  attempt  is  made  to 
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find  answers  to  some  of  the  questions  discussed  above  by 
using  a  model  which  is  one  of  the  most  advanced 
s t ress - s t r a i n  laws  presently  proposed.  The  model  was 
developed  for  sand  and  its  capabilities  were  demonstrated  by 
Lade  (1972).  Subsequent  development  by  Ozawa  (1974)  and  Wong 
(1978)  to  apply  the  model  to  geotechnical  engineering 
problems  was  not  fruitful  for  reasons  which  will  be 
discussed  as  they  arise  during  presentation. 


1 .2  OBJECTIVES 

The  research  effort  presented  here  is  directed  towards 
the  following  objectives: 

1.  To  find  out  the  significance  of  using  an  elasto-plastic 

stress-strain  law  for  modelling  the  behavior  of 

cohesionless  soils  in  a  deformation  analysis.  This  is 
accompl i shed : 

a) by  reviewing  critically  the  characteristics  and  the 
limitations  of  the  selected  model  as  well  as  the 
modifications  introduced  subsequently  by  others. 

b) by  analyzing  a  well  documented  boundary  value  problem 
using  the  finite  element  method. 

2.  To  investigate  the  reasons  why  the  previous  finite 
element  applications  of  the  model  to  geotechnical 
problems  were  not  successful. 

3.  To  compare  the  measured  collapse  load  for  a  passive 
earth  pressure  problem  against  the  predictions  of  a 
finite  element  analysis.  The  present  day  finite  element 
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models  are  not  capable  of  simulating  the  propagation  of 
a  thin  rupture  zone  in  a  previously  undetermined 
direction.  In  the  selected  model  the  soil  mass  which  has 
discrete  rupture  surfaces  at  failure  is  approximated  by 
a  continuum  with  averaged  deformation  properties. 

1  .  3  ORGANIZATION  OF  THESIS  ' 

In  Chapter  2,  the  original  formulation  of  Lade's  model 
and  the  modifications  suggested  by  others  are  introduced. 

Chapter  3  is  partly  devoted  to  the  re-evaluation  of  the 
modified  form  of  the  model.  A  shortcoming  of  the  model  is 
indicated  and  a  new  way  of  interpreting  Lade's  original 
formulation  is  proposed  as  part  of  the  present  study. 

Although  it  is  not  used  directly  in  any  of  the  analyses 
introduced  in  this  thesis,  Lade's  strain-softening  model  is 
presented  in  Chapter  4.  One  reason  for  its  presentation  is 
that,  part  of  it  is  utilized  in  Chapter  5  in  connection  with 
the  work  hardening  model.  Also  the  problems  attached  to  the 
work  hardening  laws  of  both  models  are  found  to  be  common. 
Therefore  a  clarification  is  considered  to  be  useful  for 
future  work. 

The  significance  of  a  cap-type  yield  surface  is 
elaborated  in  Chapter  5.  The  cap  yield  surface  used  by  Lade 
in  his  strain  softening  model  has  been  selected  for  use  in 
his  work  hardening  model,  which  had  only  a  cone  shaped  yield 
surface.  By  this  means  one  of  the  major  criticisms 
concerning  unsatisfactory  modelling  under  increasing 
hydrostatic  stresses  has  been  eliminated. 
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The  development  of  a  finite  element  program  using 
Lade's  work  hardening  model  without  a  cap  is  introduced  in 
Chapter  6. 

Chapter  7  deals  with  the  predictions  of  the  finite 
element  analysis  of  a  boundary  value  problem.  Results 
obtained  by  using  finite  element  method  are  compared  with 
the  passive  earth  pressure  tests  carried  out  by  Wong  (1978). 

Finally,  in  Chapter  8  the  conclusions  of  this  research 
are  stated  and  the  recommendations  for  future  research  are 
presented . 
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CHAPTER  2 


LADE'S  ELASTO-PLASTIC  WORK  HARDENING  MODEL 

2.  1  INTRODUCTION 

The  model  was  developed  by  Lade  (1972)  to  account  for 
several  aspects  of  the  stress-strai n  and  strength 
char acter i s t i cs  of  cohesionless  soils  under  general 
three-dimensional  stress  conditions.  To  a  certain  extent,  it 
has  the  capability  for  modeling: 

1 .  Non  1 i near i ty 
2  .  Shear  di 1 atancy 

3 .  Inelasticity 

4.  Stress  path  dependency 

5.  Influence  of  intermediate  principal  stress 

6.  Coincidence  of  strain  increment  and  stress  increment 
axes  at  low  stress  levels  with  transition  to  coincidence 
of  strain  increment  and  stress  axes  at  high  stress 
levels. 

The  theory  is  based  on  experimental  data  from  cubical 
triaxial  tests  on  sand,  and  uses  the  concepts  of  the 
incremental  theory  of  plasticity. 

In  the  following  sections,  the  theory  , as  introduced  by 
Lade  and  Duncan  (1975),  is  summarized,  and  the  degree  of 
success  in  its  application  is  reviewed.  The  limitations  of 
the  model  and  modifications  to  overcome  these  limitations 
are  also  indicated. 
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2.2  ELASTO-PLASTIC  STRESS-STRAIN  RELATIONSHIP 

For  the  purpose  of  the  development  of  this  theory,  it 
is  assumed  that  the  total 
divided  into  an  elastic  part 
as  shown  i n  Eqn .  2.1. 


strain  increments  id€;;  1  can  be 


d€jj  i  and  a  plastic  part  < 


Each  part  is  then  calculated  separately.  A  schematic 
illustration  of  this  concept  is  given  in  Figure  2.1  for  a 
triaxial  compresion  test. 

2.2.1  Elastic  Strain  Increments 

The  elastic  strain  increments  are  calculated  from 
Hooke's  law,  using  the  unloading-reloading  modulus,  Eur  , 
defined  by  Duncan  and  Chang  (1970): 


where  Kur  and  n  are  elastic  unloading-reloading  parameters 
(dimensionless)  and  pQ  is  the  atmospheric  pressure  expressed 
in  the  same  units  as  Eur  and  (J3  .  The  evaluation  of  these 
parameters  from  triaxial  tests  is  explained  in  Appendix  A. 

The  value  of  Poisson's  ratio  is  assumed  to  be  zero. 
This  assumption  is  based  on  triaxial  test  data  where  the 
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axial  and  volumetric  strains  are  very  nearly  equal  for  the 
first  increment  of  load. 

As  a  result,  the  incremental  stress-s trai n  relationship 
for  an  elastic  isotropic  material  takes  the  following  form: 
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Eqn .  2.3 


2.2.2  Plastic  Strain  Increments 

For  the  calculation  of  plastic  strain  increments,  the 
model  follows  the  basic  requirements  of  plasticity  theory  as 
summarized  by  Lade  and  Duncan  (1975)  : 

"There  are  three  basic  requirements  for  a  plastic 
stress-strai n  theory: 


1.  There  must  exist  a  yield  surface  such  that  if 
the  soil  is  subjected  to  changes  in  stress 
represented  by  points  inside  that  surface,  the 
soil  will  deform  elastically,  whereas  if  the 
changes  in  stress  tend  to  cross  the  yield 
surface,  it  will  simultaneously  yield 

plastically  and  deform  elastically.  The  yield 
surface  expands  as  the  soil  is  loaded  to 


< 

1 

, 


successively  higher  stress  levels,  and  at 
failure  the  yield  surface  coincides  with  the 
f ai lure  surface . 

2.  A  flow  rule  is  also  required.  The  flow  rule  is  a 

law  which  relates  the  relative  magnitudes  of  the 
strain  increments  to  the  stresses.  The  flow  rule 
is  derived  from  the  requirement  that  the  plastic 
strain  increment  direction  should  be  normal  to 
the  plastic  potential  surface.  There  are  an 
infinite  number  of  such  surfaces,  and  one  passes 
through  every  point  in  stress  space.  In 

so-called  "perfect  plasticity,"  the  plastic 
potential  surface  is  assumed  to  be  the  same  as 
the  yield  surface.  As  mentioned  previously,  this 
implies  much  higher  rates  of  dilation  than  are 
observed  in  tests  on  real  soils.  When  the 
plastic  potential  surface  is  assumed  to  differ 
from  the  yield  surface,  i.e.  when  the  flow  rule 
is  "non-associ ated , "  much  better  agreement 
between  theory  and  experiment  can  be  achieved. 

3.  A  work-hardeni ng  law  is  needed,  from  which  the 
magnitudes  of  the  plastic  strain  increments 
caused  by  a  given  stress  increment  can  be 
determi ned . " 
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The  details  of  Lade's  model  in  relation  to  those 
requirements  are  given  below: 

1 .  Yield  surface 

The  yield  surface  is  assumed  to  have  the  shape  of  a 
cone  with  the  apex  at  the  coordinate  center  of  the 
principal  stress  space  as  shown  in  Figure  2.2.  It  is 
expressed  as  a  function  of  the  first  and  third  stress 
invariants  as  follows: 

T3 

f  =  -±L  .  Eqn.  2.4 

i3 

f  denotes  the  stress  level  and  its  value  varies  from  27 
for  hydrostatic  stress  conditions  up  to  a  value  at 
failure.  With  increasing  values  of  f,  the  yield  surface 
expands  continuously,  and  becomes  identical  with  the 
failure  surface  at  its  outermost  shape. 

2.  Flow  Rule 

A  non-associ a  ted  flow  rule  is  used  for  this  model 
in  which  the  plastic  potential  surface  is  no  longer 
assumed  to  be  the  same  as  the  yield  surface.  In  doing 
so,  much  better  predictions  for  the  rate  of  dilation  are 
expected . 

The  plastic  potential  is  a  function  of  the  state  of 
stress,  from  which  the  relative,  but  not  the  absolute 
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magnitudes  of  the  components  of  the  plastic  strain 
increments  can  be  determined  by  differentiating  with 
respect  to  stresses.  The  function  incorporated  in  the 
theory  is  expressed  by  Lade  and  Duncan  (1975)  in  a  form: 


9  =  I?  -  K2  *  I3 


Eqn.  2.5 


where : 

g  is  the  value  of  the  plastic  potential  and 
K2  is  a  constant  for  a  given  value  of  f.  It  is 
calculated  from 


K2  =  Axf  ♦27x(1-A) 


Eqn.  2.6 


A  is  a  material  constant  and  its  evaluation  from 
triaxial  compression  tests  is  given  in  Appendix  A. 

Eqn.  2.5  describes  a  series  of  surfaces  which  are  normal  to 
the  plastic  strain  increment  directions. 

In  the  theory  of  plasticity,  see  for  example  Hill 
(1950),  the  relation  between  the  plastic  strain  increments 
and  the  plastic  potential  function  is  given  by  an  expression 
as  fol lows : 


Eqn.  2.7 


in  which  X  is  called  the  proportionality  constant.  The 

is  the  same  for  all  components  of  plastic 


value  of 
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strain  increments  and  is  determined  by  the  work  hardening 
law  as  explained  next. 


3.  Work  hardening  law 

The  work  hardening  law  adapted  for  the  model  is  an 
experimentally  determined  relationship  between  the  plastic 
work  Wp  and  the  stress  level  f.  The  increment  of  plastic 


work  done  per  unit  volume  due  to  a  strain  increment  < 

|T 


dSf; 


G-;  <jd  €|j  \  and  is  added  to  the  previous 


is  calculated  from  < 
values  to  find  the  total  plastic  work. 

Test  data  by  Lade  (1972)  have  indicated  that  the 
plastic  work  was  very  small  for  a  range  of  f  starting  with 
27  at  hydrostatic  state  of  stress  up  to  a  certain  value 
which  was  called  the  threshold  stress  level  f^  .  As  a 
convenience  in  fitting  curves  to  the  experimental  data,  it 
has  been  assumed  that,  for  valus  of  f  between  27  and  f^  ,  no 
plastic  strains  occur  and  no  plastic  work  is  done. 
Therefore,  the  relations  between  Wp  and  (f-f^  )  have  been 
approximated  by  hyperbolae  for  which  Eqn.  2.8  has  been 
proposed . 


Wp 

a  +  b  x  Wp 


Eqn.  2.8 


The  initial  slope  of  a  curve  representing  Wp  versus  (f-ft  ) 
relationship  is  the  reciprocal  of  the  parameter  a  as  shown 
in  F i gure  2.3. 

The  value  of  a  increases  with  confining  pressure, 


and  this 


' 
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variation  has  been  expressed  as  in  Eqn.  2.9. 


a  =  M  x 


Eqn .  2.9 


where : 

p  is  atmospheric  pressure  express  the  same 

a 

units  as  a  and  (J^  . 

M  and  1  are  dimensionless  numbers.  Their  evaluation 

from  test  data  is  given  in  Appendix  A. 

The  parameter  b  in  Eqn  2.8  is  the  reciprocal  of  the  ultimate 
value  of  (f-fj  )  which  the  hyperbola  approaches 
asymptotically  with  increasing  values  of  Wp  .  This 
relationship  is  given  by  Eqn .  2.10. 


1  .  Eqn.  2.10 


Since  the  value  ( f - f t  )u^  determined  from  the  curve  fitting 
procedure  is  always  larger  than  the  value  of  (f-f^  )  at 
failure  for  all  finite  values  of  Wp ,  a  new  parameter  rf  was 
introduced  to  relate  the  asymptotic  value  of  (f-ft  )  to  its 
value  at  fai lure. 

(K^ 
f  (f-fdult 


Eqn.  2.11 
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The  determination  of  the  proportionality  constant  A 
follows  the  development  outlined  by  Hill  (1950)  and  is  given 
as : 


x  =  d  Wp 

39 


Eqn .  2.12 


The  increment  of  plastic  work  dWp  is  expressed  as: 


axdf 


(i-rfxizfLY 


Eqn.  2.13 


The  plastic  strain  increments  expressed  in  suffix 
notation  in  Eqn.  2.7,  can  now  be  written  in  matrix  form: 


'  d€$  ' 

3if-K2(cryaz-Ty2z) 

d  €y 

3  II  “  K2(  0-zCJx- Tzx ) 

d€z 

✓ 

axdf 

3ll  “  K2(0^  CJy  -  7^y  )  ^ 

d€£y 

n 

2K2(azTxy-TyzTzx) 

d  €yz 

2  K2(axTyz-TxyT2X) 

d£pz* 

* 

2  K2((TyTzx- TyZ  7^y) 

Eqn.  2.14 


Therefore,  the  elastic  and  plastic  parts  of  strain 
increments  are  fully  described.  Substitution  of  Eqn.  2.3  and 


' 
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Eqn .  2.14  into  Eqn .  2.1  allows  the  total  strain  increments 
to  be  calculated  if  the  state  of  stress  and  the  stress 
increments  are  Known. 

There  are  nine  soil  parameters  required  to  characterize 
the  soil  behavior  according  to  this  theory,  and  their  values 
can  be  derived  from  the  results  of  conventional  triaxial 
compression  tests. 

2.3  LIMITATIONS 

The  theory,  as  recapitulated  in  the  previous  section, 
is  not  able  to  model  all  aspects  of  soil  behavior. 

1.  Strain  softening  cannot  be  handled.  If  the  conditions 
for  strain  softening  specified  in  Chapter  4  prevail, 
then  the  predictions  may  deviate  from  measured  results 
considerably  depending  on  the  degree  of  softening. 

2.  According  to  the  theory,  no  plastic  straining  takes 
place  as  a  result  of  a  proportional  loading.  This  is  due 
to  the  fact  that  the  model  does  not  have  a  yield  cap. 
For  example  in  hydrostatic  compression  tests  which 
follow  a  proportional  loading  path,  the  model  predicts 
recoverable  strains  only.  This  is  in  contradiction  with 
measured  behavior.  The  degree  of  discrepancy  between  the 
predicted  and  measured  results  depends  mainly  on  the 
relative  density  of  the  soil  and  the  stress  level. 

3.  Although  the  Mohr  envelope  is  curved  for  most 

cohesionless  soils,  the  failure  envelope  in  the  triaxial 
plane  for  this  model  is  linear.  The  effects  of  this 


. 
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limitation  apply  more  to  soils  under  low  mean  normal 
stresses  where  the  curvature  of  failure  envelope  is 
appreciable . 

4.  The  value  of  Poisson's  ratio  for  elastic  deformations 
was  assumed  to  be  equal  to  zero.  Therefore  the  terms 
multiplied  by  its  value  will  disappear  in  formulae  to 
calculate  strain  increments.  Consequently,  no  straining 
should  be  expected  in  any  other  direction  but  that  of 
the  stress  increment. 

2.4  APPLICATIONS  OF  THE  MODEL 

Applications  of  the  model  to  laboratory  tests  and 

engineering  problems  are  reviewed  next. 

1.  Predictions  of  laboratory  testing  results: 

To  assess  the  ability  of  the  model  to  predict  soil 
behavior,  a  variety  of  tests  with  different  stress  paths 
were  conducted  and  the  comparative  results  were 
published  by  Lade  and  Duncan  (1975,1976).  Their 
conclusions  were: 

"The  generally  good  agreement  between  the 
calculated  and  measured  strains  indicated  that 
the  conditions  of  primary  loading,  unloading, 
and  reloading  implied  in  plasticity  theory  for 
isotropically  worK-hardeni ng  materials  is 
reasonably  accurate  for  cohesionless  soils.  The 
strains  calcul ated . . . are  in  reasonable  agreement 
with  those  measured,  but  some  discrepancies  do 


' 
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occur.  The  largest  differences  between  measured 
and  calculated  strains  occur  for  proportional 
loading  with  increasing  stresses  and  for 
unloading  and  reloading  at  constant  confining 
pressure.  Considering  that  the  strains  observed 
in  these  cases  are  relatively  small,  the 
calculated  strains  may  be  close  enough  for  many 
purposes.  The  analyses  thus  demonstrate  the 
usefulness  of  el astopl ast ic  stress  strain  theory 
for  calculating  strains...  during  complex  stress 
changes  ..." 

The  effectiveness  of  the  theory  in  modeling  stress 
strain  behavior  along  several  different  stress  paths  was 
also  confirmed  by  Medeiros  (1979). 

Engineering  Problems: 

Ozawa  (1973)  was  the  first  to  use  the  model  for  the 
analysis  of  stresses  and  movements  in  earth  masses.  His 
results  of  finite  element  analyses  of  passive  earth 
pressure  problems  were  published  later  by  Ozawa  and 
Duncan  (1976a).  Due  to  errors  in  the  analysis  of  this 
published  data,  and  the  problems  related  to  the 
formulation  of  the  model,  (see  Section  2.5),  the  results 
were  not  suitable  for  determining  the  usefulness  of  the 
mode  1 . 

The  second  attempt  to  find  the  effectiveness  of  the 
model  for  engineering  problems  was  made  by  Wong  (1978). 
Using  the  modifications  summarized  in  the  following 
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section,  the  predictions  of  a  finite  element  analysis 
were  compared  with  experimental  results  from  passive 
pressure  tests  on  sand.  The  calculated  values  of  average 
stresses  on  the  wall  exceeded  those  measured  by  a 
considerable  amount.  Attention  was  then  drawn  to  the 
limitations  of  the  model  as  possible  reasons  for  the 
discrepancies.  For  the  same  problem,  Wong  has  also 
provided  calculated  results  obtained  by  using  simpler 
relationships,  such  as  the  nonlinear  elastic  hyperbolic 
model.  The  predictions  made  by  using  Lade's  model  were 
not  as  impressive  as  the  others,  yet  the  simpler 
stress-strai n  relationships  have  more  severe 
1 imi tat  ions . 


2.5  PREVIOUS  MODIFICATIONS 

Modifications  and  corrections  for  the  model  were  first 
proposed  by  Ozawa  and  Duncan  (1976b)  as  follows: 


1.  The  plastic  potential  function 

was  found  to 

be 

mathemat i 

cally  inconsistent  with 

the 

express i ons 

for 

plastic 

strain  increments  and 

a 

new  equation 

was 

proposed 

as : 

g  =  I3x(f-27 

.  Eqn . 

2.  15 

where : 


' 
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OL  is  a  constant  for  a  particular  material  and  is 
defined  as: 


where  A  is  the  same  as  in  Eqn .  2.6 

2.  Due  to  the  error  in  the  calculation  of  derivatives  of 
plastic  potential  function  with  respect  to  shear 
stresses,  the  shear  strain  increments  were  in  an  error 
by  a  factor  of  2.0.  Therefore  on  the  right  hand  side  of 
Eqn.  2.14,  the  multiplier  2.0  of  each  term  related  to 
plastic  shear  strain  increments  has  to  be  removed. 

3.  The  parameter  a  of  Eqn.  2.8  varies  with  (J '  according  to 
Eqn.  2.9.  This  variation  was  neglected  in  calculating 


df/d  O'  which  is  used  in  the  finite  element  formulation 
to  form  the  elasto-plast ic  constitutive  matrix 


Having  made  these  corrections,  the  authors  then 
presented  the  elasto-plastic  constitutive  relation  in 
matrix  notation,  suitable  for  use  in  a  finite  element 
analysis.  Wong's  predictions  were  based  on  this  form  of 
the  model  and  therefore  utilized  all  previous 
modi f icat ions . 


The  same  corrections  were  again  published  by  Duncan 
et  al  (1977)  with  an  illustration  of  the  new  plastic 
potential  surface  which  has  an  unsymmetr i ca 1  bullet 
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shape  in  principal  stress  space.  The  intersections  of 
the  plastic  potential  surfaces  with  the  triaxial  plane 
are  shown  in  Figure  2.4. 


27 


Figure  2.1  Elastic  and  plastic  parts  of  strain  in  triaxial 
test  (After  Lade  and  Duncan  (1975)) 
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Figure  2.2  Yield  and  failure  surfaces  in  (a)  principal 
stress  space  (b)loci  of  surfaces  on  octahedral  plane. 
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Figure  2.3  Relationship  between  plastic  work  and 
level.  (After  Ozawa  and  Duncan,  1976a) 
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Figure  2.4  Trace  of  plastic  potential  surfaces  on  tr 
plane  (After  Duncan  et  a  1 , 1977) 
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CHAPTER  3 


RE-EVALUATION  OF  WORK  HARDENING  MODEL 
3. 1  INTRODUCTION 

For  most  stress  paths,  Lade's  model  has  been  shown  to 
be  effective  in  its  predictions  of  stress  strain  behavior  of 
cohesionless  soils  in  laboratory  tests.  However  finite 
element  applications  of  the  model  have  not  reached  the  same 
level  of  success  so  far.  In  the  past,  difficulties  were 
encountered  with  regard  to  the  formulation  of  the  model,  and 
consequently  a  number  of  modifications  were  proposed  by 
others  as  reviewed  in  Section  2.5. 

The  discouraging  results  of  finite  element  analyses  of 
earth  pressure  problems  by  Wong  (1978)  as  opposed  to  the 
satisfactory  predictions  of  laboratory  tests  by  Lade  (1972), 
Lade  and  Duncan  (1976),  suggested  that  either  the 
limitations  of  the  model  were  too  severe  or  the 
modifications  proposed  by  Ozawa  and  Duncan  (1976b)  have  not 
been  effective  enough  for  the  model  to  become  useful  in 
engineering  practice. 

In  this  Chapter,  certain  aspects  of  the  original 
formulation  of  the  model  and  previous  modifications  are 
examined.  The  significance  of  correct  interpretation  of  the 
laws  used  in  Lade's  original  theory  is  demonstrated  by 
comparing  the  test  results  of  Lade  (1972)  with  predictions 
by  the  previously  modified  and  presently  re-evaluated  forms 
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of  the  model . 

The  adaptation  of  the  suggestions  made  here  to  finite 
element  analysis  is  dealt  with  in  Chapter  6. 

3.2  VALIDITY  OF  PREVIOUS  MODIFICATIONS 

Any  alteration  to  the  original  formulation  of  the  model 
affects  its  ability  to  predict  stress-strain  behavior  of 
soils.  To  evaluate  the  effect  of  modifications  on  the 
existing  capabilities  of  the  model,  predictions  have  to  be 
compared  with  laboratory  test  results.  The  calculations  made 
by  using  the  original  form  of  the  model  (Lade,  1972)  and  the 
form  with  previous  modifications  (Ozawa  and  Duncan,  1976) 
are  presented  next  with  some  test  data. 

Ozawa  and  Duncan  (1976)  have  proposed  that  the 
variation  of  parameter  a  with  respect  to  0"3  has  to  be 


used  in 


included  in  the  calculation  of 


the  finite  element  formulation  to  form  the  elasto-plast ic 


If  this  proposal  is  used  in 


constitutive  matrix 


the  finite  element  analysis,  the  results  will  be  equivalent 
to  that  of  a  simple  stress  analysis  where  the  parameter  a  is 
used  as  a  variable  in  the  calculation  of  plastic  work 
increments.  Therefore,  realizing  that  the  parameter  a  of 
Eqn .  2.8  has  to  be  treated  as  a  variable  rather  than  a 
constant  and  using  the  new  plastic  potential  function  given 
by  Eqn.  2.15,  the  constitutive  relationship  in  modified  form 


can  be  wr i t  ten  as : 


. 


Eqn .  3.1 


d€xp 

3i?-x2((jy<7z  _ry2z) 

d  €y 

3I?-K2(crz(7x  -Tzl  ) 

d€j? 

3i?-'<2(axay-Tx2y) 

d€py 

>  =  MULT x  < 

K2(Txyaz-TyzTzx)  r 

d  £yZ 

k2(7^zO"x  -  Tzx  Txy  ) 

d€zx 

. 

K2(Tzxay-  TxyTyz) 

* 

where : 


MULT 


(M*bhftcr,h  d  fl+MicrjVdcr,) 


3Mg 


Q!-1 

la(f-27) 


k2=  a,  f  +  a2 
a,  =  i  -  MOL 
A2=  21/01 
Oi-  1/(1  -  A) 
h  =  W  p  /  ( 0"3  ) 

Notes:  1)  Wp  has  to  be  updated  at  each  loading  step 

2)  The  column  vector  on  the  right  hand  side  of  Eqn. 
3.1  will  be  denoted  by  |cOLUMn|  and  the  left  hand 
side  by|d€^|in  the  following  as  required. 

3)  M,  g,  b,  1  are  all  defined  in  Chapter  2. 

Several  stress-strain  curves  following  different  stress 
paths  were  reproduced  by  substituting  Eqn.  3.1  into  Eqn. 
2.1,  and  are  presented  with  the  test  results,  if  available, 
in  Figures  3.1  to  3.5.  The  predictions  by  using  Lade's 
approach  are  also  plotted  in  the  same  figures.  Although  the 
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predicted  curves  of  the  previously  modified  approach  are  as 
good  as  the  curves  produced  by  using  Lade's  approach  for 
Figure  3.1  and  3.2,  the  gap  between  these  curves  rapidly 
opens  up  for  large  changes  in  (J3  but  relatively  small 
changes  in  f.  Figures  3.3  to  3.5  demonstrate  that  the 
modified  approach  may  produce  significantly  different 
results  from  Lade's  approach,  even  resulting  in  different 
signs  for  strains  as  shown  in  Figure  3.5. 

Although  there  are  no  measured  test  results  to  be 
plotted  in  Figures  3.4  and  3.5  to  support  either  one  of  the 
predictions,  it  can  be  concluded  that  the  use  of  Eqn.  3.1 
alters  the  abilities  of  the  model  for  certain  stress  paths. 

Figure  3.6  shows  the  strain  increment  vectors  for  two 
different  stress  paths  in  the  triaxial  plane  as  predicted  by 
the  previously  modified  form  and  the  presently  re-evaluated 
approach  which  is  described  later  in  this  chapter. 


For 

s i mi  1 ar 

stress  paths, 

Lade  and 

Duncan  (  1976) 

provided 

measured 

test  results 

which  are 

reproduced  on 

Figure  3.7.  A  comparison  of  measured  and  predicted  strain 
increment  vectors  shows  that  the  re-evaluated  procedure 
predicts  the  increments  of  strain  reasonably  well  for  both 
stress  paths,  while  the  modified  approach  fails  to  do  so  for 
the  stress  path  with  decreasing  deviator  stress  and 
decreasing  confining  pressure. 

To  find  out  which  one  of  the  previous  modifications  is 
responsible  for  altering  the  abilities  of  the  model,  the 
effect  of  each  modification  to  the  original  formulation  is 
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investigated  separately.  The  modification  related  to  the 
plastic  potential  function  is  treated  first,  with  the 
parameter  a  held  constant  in  calculating  the  plastic  work 
increment.  Rather  than  reproducing  the  s tress-strai n  curves, 
the  following  proof  is  provided  to  show  that  the  modified 
plastic  potential  function  (Eqn.  2.15)  does  not  change  the 
stress-strain  relation  in  incremental  form. 

To  differentiate  between  the  original  and  modified 
plastic  potential  functions,  g  of  Eqn.  2.5  and  Eqn.  2.15  are 
renamed  here  as  g0  and  gm  respectively.  Substituting  Eqn. 
2.12  into  Eqn.  2.7  and  using  Eqn.  2.5  for  the  plastic 
potential  function  results  in  the  following  plastic  strain 
increment  values: 


COLUMN 


Eqn .  3.2 


I 


where  :<jd€  ?and  jCOLUMN 


*  are  defined  in  the  notes  for  Eqn 


3.1. 

If  Eqn.  2.15  is  used  rather  than  Eqn.  2.5,  the  set  of 
equations  for  calculating  the  plastic  strain  increments 
becomes : 


d€p 


V39m 


< COLUMN  } 


Eqn .  3.3 


If  the  multipliers  on  the  righthand  side  of  Eqn.  3.2  and 
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Eqn.  3.3  have  the  same  values,  then  d  £P.  calculated  by 

|J 

either  way  are  identical.  This  can  readily  be  shown: 

Ignoring  the  common  multiplier  (dWp/3.0),  and  writing 
f,  gQ  and  gm  in  terms  of  stress  invariants  gives: 

Reduced  multiplier  for  Eqn.  3.2: 


1 


I?-K2I3 


1 


if-tA  il  +  27  (1  -  A)]  I3 
I3 


J_  _  _ ] _ 

9°  if- A  1^-27 13  +27AI3 


Eqn.  3.4 


Reduced  multiplier  for  Eqn.  3.3: 


CM 


9m 


Ot  (f-  27) 


OM 

1  {f- 27) 


1 


1  -A  I3(  f  —  2  7) 


a 


( 1  -  A)  (  f  -27)  I. 


9m 


at- 1 

at  (f- n) 


1 


if  -  Alf  -  2  7 13  +27A I3 


Eqn.  3.5 


where  O'  =  1/(1  -  A ) 

Eqn.  3.4  and  Eqn.  3.5  have  equal  values,  therefore,  the 
use  of  Eqn.  2.15  for  the  modified  plastic  potential  function 
will  produce  results  identical  to  those  obtained  from  Lade's 
original  formulation.  If  this  modification  is  used,  it 
improves  the  formulation  of  the  model  in  the  sense  of 
mathematical  integrity.  The  explanation  for  this  statement 
is  as  follows:  when  the  derivatives  of  the  plastic  potential 
function  are  needed  for  the  flow  rule,  the  use  of  Eqn.  2.5 
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requires  that  K2  be  held  constant  in  the  derivation 
process.  Yet  K2  is  a  function  of  f  which  is  in  turn  a 
function  of  stresses.  The  use  of  Eqn.  2.15  does  not  require 
such  an  assumption. 

( 

In  the  light  of  the  above  observations  regarding  the 
plastic  potential  function,  the  cause  for  large 
discrepancies  between  the  predictions  by  Lade's  approach  and 
the  results  of  the  modified  formulation,  as  plotted  in 
Figure  3.3  to  Figure  3.6,  seems  to  be  related  to  the 
parameter  a,  in  other  words,  to  the  work  hardening  law  of 
the  theory. 


3.3  PROPOSED  INTERPRETATION  OF  WORK  HARDENING  LAW 


The  work  hardening  law  and  its  modification  are 
re-evaluated  here  to  establish  a  correct  method  for  its  use. 

Lade's  test  results  from  which  the  work  hardening  law 
was  obtained  are  reproduced  in  Figure  3.8.  The  results  show 
that  there  is  a  unique  relationship  between  f  and  Wp  for 
tests  where  0~3  is  held  constant.  This  relationship  is 
expressed  by  Lade  and  Duncan  (1975)  in  mathematical  form 
(see  Eqn.  2.8).  Figure  3.8  also  shows  some  data  points  for 

C72  -  (J3  ^ 

various  b  values  [b  =— — —  ]  obtained  from  tests  where  U3 


was  still  kept  constant.  These  results  allowed  Lade  and 
Duncan  (1975)  to  state  that  "the  relationship  between  Wp  and 
f  ...  depends  on  the  confining  pressure,  (J3  ,  but  it  is 
essentially  the  same  for  all  values  of  b" .  There  are  no 
other  test  results  with  plastic  work  calculations  to  support 


. 
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the  validity  of  Eqn .  2.8  for  stress  paths  with  varying  (J3  • 

Therefore,  this  equation  is  applicable  strictly  for  stress 
changes  for  which  CJ3  remains  constant. 

In  order  to  use  Eqn.  2.8  in  the  analysis  of  a  general 
three  dimensional  problem  where  ( J3  generally  changes 
during  loading,  an  assumption  has  to  be  made.  The  theory 
which  was  first  published  by  Lade  and  Duncan  (1975)  does  not 
elaborate  on  this  point,  but  simply  treats  the  parameter  a 
as  a  constant  when  Eqn.  2.8  is  used  to  derive  an  expression 
for  the  increment  of  plastic  work.  However,  in  an  attempt  to 
use  Eqn.  2.8  in  its  true  sense,  Lade  (1972)  calculated  the 
stress  dependent  parameter  a  for  an  average  CJ3  value 
which  was  obtained  from  initial  and  final  values  of  (J3 
for  a  given  stress  increment.  Some  of  his  predictions  of 
laboratory  tests  can  be  found  in  the  publication  by  Lade  and 
Duncan  (1976).  Unfortunately  the  importance  of  the 

non-generality  of  Eqn.  2.8  for  changing  C73  values  has 
never  been  emphasized,  nor  has  the  basic  assumption  which 
allows  its  use  in  such  problems  been  clearly  stated.  In  the 
publications  by  Ozawa  and  Duncan  (1976,  1976b),  Eqn.  2.8 

remained  as  a  three  dimensional  representat ion  of  the  work 
hardening  law.  Subsequently,  as  a  deviation  from  Lade's 
original  work,  Duncan  et  al .  (1977)  and  Wong  (1978)  treated 
the  parameter  a  as  a  variable  when  a  differentiation  was 
needed  on  Eqn .  2.8. 

The  consequences  of  assuming  that  Eqn.  2.8  is  unique 
for  general  three  dimensional  problems  are  as  follows: 
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Case  l: 

Referring  to  Figure  3.9a,  the  initial  state  of  stress 
is  assumed  to  be  on  the  current  yield  surface  as  indicated 
by  point  A  with  a  total  plastic  work  value  equal  to  WA .  An 
increment  of  stress  is  applied  such  that  the  stress  path 
follows  the  route  from  A  to  B  without  travelling  outside  the 
current  yield  surface.  The  limits  of  the  elastic  region  are 
not  exceeded,  so  there  are  no  plastic  deformations,  and  no 
increment  in  plastic  work.  Therefore,  the  total  plastic  work 
done  is  the  same  for  points  A  and  B.  But  since  the  value  of 
CJ3  has  changed,  and  with  it,  the  value  of  parameter  a,  the 
Eqn .  2.8  is  no  longer  satisfied  for  point  B. 

Case  1 1 : 

Figure  3.9b  illustrates  three  different  stress  paths, 
all  of  which  start  from  the  hydrostatic  axis.  All  three 
follow  the  stress  path  of  the  conventional  triaxial  test  up 
to  the  new  yield  surface,  then,  without  leaving  that  surface 
end  at  point  F.  The  change  in  f  value  in  going  from  the 
hydrostatic  axis  to  the  new  yield  surface  is  the  same  for 
all  points  A,  B  and  C.  But  the  parameter  a  calculated  for 
point  A  wi 1 1  be  the  smallest  of  all,  and  the  plastic  work 
increment  for  path  A  to  D  wi 1 1  be  smaller  than  for  path  B  to 
E  or  C  to  F .  Therefore,  depending  on  which  path  is  followed, 
the  total  plastic  work  at  point  F  wi 1 1  have  a  different 
value.  Although  all  these  work  values  as  well  as  an  infinite 
number  of  others  are  possible  for  point  F,  only  one  of  them 
satisfies  Eqn.  2.8.  Subsequent  increments  of  stress  from 
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point  F  with  increasing  stress  level  will  create  problems 
for  any  total  plastic  work  value  other  than  the  one  which 
satisfies  Eqn.  2.8.  The  reason  for  this  is  that  Eqn .  2.8 
will  become  an  inequality  and  cannot  be  used  any  further  to 
calculate  the  increments  of  plastic  work,  which  are  needed 
for  Eqn.  2.12. 

Case  III: 

This  case  is  related  to  the  use  of  Eqn.  2.8  in  its 
modified  form,  in  which  the  parameter  a  is  treated  as  a 
function  of  (J3  in  calculating  the  increment  of  plastic 
work,  as  formulated  below: 


d  WD=— 
p  ML 


(M+bhMtjJ  )( df  )+Mloi"'h(dCr,) 


. . .  Eqn .  3.6 


The  problem  does  not  surface  with  this  formulation  until  the 
term  with  d  (J3  becomes  larger  in  absolute  value  with  a 
negative  sign  than  the  term  with  df  on  the  right  hand  side 
of  the  equation  as  in  stress  paths  with  decreasing  (J3  and 
increasing  f  values.  As  a  result  of  this,  d  Wp  will  have  a 
negative  value.  As  long  as  the  plastic  deformation  is  an 
irreversible  process  from  which  the  energy  cannot  be 
recovered,  a  negative  d Wp  will  contradict  the  laws  of  the 
theory  of  plasticity. 

Due  to  the  lack  of  generality  as  demonstrated  by  Case  I 
and  II,  Eq.  2.8  is  not  assumed  here  as  a  work  hardening  law 
applicable  for  all  problems.  Instead  the  following  procedure 
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is  adopted: 

Rather  than  searching  for  a  general  relationship 
between  f  and  Wp  ,  and  then  differentiating  that  function  to 
find  a  relation  between  df  and  dWp  ,  the  following 
assumption  is  made  to  establish  a  criterion  about  how  the 
yield  surface  will  expand  as  the  plastic  work  is  done.  Using 
Figure  3.8,  the  tangent  to  the  related  experimental  curve 
specified  by  (J3  ,  at  the  point  correspondi ng  to  known  f, 
defines  the  relation  between  the  increments  of  f  and  Wp 
without  any  regard  for  the  actual  value  of  Wp.  As  it  is 
shown  earlier,  depending  on  the  loading  history  of  the  soil, 
the  accumulated  value  of  Wp  may  vary  for  the  same  values  of 
f  and  0"3  .  However,  the  increments  of  plastic  work  will 
always  be  unique  for  a  given  state  of  stress  and  stress 
increment.  In  mathematical  form  this  statement  is  equivalent 
to  the  expression  given  by  Eqn .  2.13,  which  is  Lade's 
original  equation  for  plastic  work  increments.  Therefore,  it 
can  be  concluded  that  the  re-evaluated  approach  makes  full 
use  of  the  experimental  data  and  Eqn.  2.8  without 
influencing  adversely  what  has  been  formulated  so  far. 

To  find  out  the  results  of  such  an  approach  the 
stress-strai n  curves  of  Figures  3.1  to  3.6  were  reproduced 
and  plotted  as  "re-evaluated"  predictions  on  the  same 
figures.  These  results  show  that  the  curves  obtained  by  the 
re-evaluated  approach  follow  closely  the  original 
predictions  which  were  provided  by  Lade  (1972)  for  the 
verification  of  his  model. 
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3.4  CONCLUSIONS 

As  the  calculations  provided  here  indicate,  the 
measured  stress-strai n  behavior  of  cohesionless  soils  can  be 
successfully  predicted  by  using  the  newly  proposed 
interpretation  without  making  the  previous  modifications  of 
Ozawa  and  Duncan  (1976).  At  the  same  time,  the  problems  that 
are  created  in  finite  element  analysis  by  using  Eqn .  2.8  as 
the  work  hardening  law  can  be  avoided  . 
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o  o  o  L»de's  approach  (prediction) 

□  □  □  Previous  modifications  (prediction) 

-  Re-evaluated  ( prediction) 

AAA  Experimental 


6,(%) 


Figure  3.1  Stress-path  with  increasing  deviator  stress 
constant  confining  pressure . Monterey  No.O  Sand 
Experimental  results  after  Lade  and  Duncan  (1976) 
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o  o  o  Lade's  approach  (prediction) 

□  □  □  Previous  modification  ( prediction ) 
-  Pe-evaluated  ( prediction) 


AAA  Experimental 


Figure  3.2  Stress-path  with  increasing  deviator  stress  and 
decreasing  confining  pressure.  Monterey  No.O  Sand. 
Experimental  results  after  Lade  and  Duncan  (1976) 
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o  o  o  Lade's  approach  (  prediction  ) 

□  □  □  Previous  modifications  (  prediction  ) 
-  Re-evaluated  (  prediction ) 


Figure  3.3  Stress-path  with  decreasing  deviator  stress  and 
decreasing  confining  pressure.  Monterey  No.O  Sand. 
Experimental  results  after  Lade  and  Duncan  (1976) 
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o  o  o  Lade's  approach  (prediction  ) 

□  □  □  Previous  modifications  ( prediction  ) 
-  Re-evaluated  (  prediction  ) 


Figure  3.4  Stress-path  with  increasing  deviator  stress  and 
increasing  confining  pressure.  Monterey  No.O  Sand. 
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o  o  o  Lade's  approach  (prediction  ) 

□  □  □  Previous  modifications  ( prediction ) 
-  Re-evaluated  (prediction) 


Figure  3.5  Stress-path  with  decreasing  deviator  stress  and 
decreasing  confining  pressure.  Monterey  No.O  Sand. 
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Figure  3.6  Predicted  strain  increment  vectors 
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VT.  a3  (kg/cm2)  ,  VT.  €3 


Figure  3.7  Measured  Directions  of  Strain  Increment  Vectors 
for  Different  Stress-Paths  in  Triaxial  Plane  (After  Lade  and 
Duncan,  1976) 
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Figure  3.8  Variation  of  Plastic  Work  for  Dense  Monterey  No.O 
Sand  (After  Lade  and  Duncan,  1976) 
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Figure  3.9  Stress  increments  with  a)  no  change  in  stress 
level  b)  conventional  triaxial  test  path  followed  by  no 
change  in  stress  level. 


CHAPTER  4 


STRAIN  SOFTENING  MODEL 


4.  1  INTRODUCTION 

In  his  continuing  work  on  soil  plasticity  Lade  (1977) 
developed  a  strain  softening  model  to  overcome  the 
limitations  of  his  work  hardening  model.  Those  limitations 
are  presented  in  Section  2.3.  Although  some  of  the  general 
concepts  used  for  the  work  hardening  model  remained  the 
same,  the  formulation  of  the  conical  yield  surface,  the 
corresponding  plastic  potential,  and  the  hardening  law  were 
all  modified  in  the  development  of  the  strain  softening 
model . 

In  this  model  the  cone  has  a  curvature  like  a  bullet 
and  not  only  expands  as  work  hardening  takes  place  but  can 
also  contract  in  order  to  model  the  strain  softening  after 
the  failure  surface  is  reached.  Since  the  failure  surface  is 
the  outermost  yield  surface  which  has  a  curvature,  the 
failure  envelope  in  the  triaxial  plane  is  no  longer  straight 
but  curvi 1 i near . 

In  addition  to  the  strain  softening  feature,  a  cap  type 
yield  surface  was  added  to  the  open  end  of  the  conical  yield 
surface  to  account  for  the  plastic  strain  increments  due  to 
proportional  loading. 

Elastic  strains  were  calculated  using  a  non-zero 
Poisson's  ratio  (zero  was  used  for  the  work  hardening  model) 
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and  therefore  a  general  improvement  was  expected  in 
predictions  of  elastic  deformations. 

The  strain  softening  model  has  been  used  so  far  mainly 
for  the  analysis  of  laboratory  tests.  The  modelling  of 
stress-strain  curves  was  reasonably  accurate  and  the  pore 
pressure  calculations  in  undrained  loading  were 
satisfactory.  The  model  was  developed  for  sand  but  its 
application  to  cohesive  soils  was  also  a  success  (Lade  and 
Musante,  1976).  The  finite  element  application  of  this  model 
has  not  been  popular  yet  even  among  researchers.  At  present, 
the  only  published  data  is  by  Aubry  and  Des  Croix,  (1979), 
and  the  correctness  of  their  results  is  not  beyond  doubt. 

The  present  research  work  essentially  deals  with  the 
work  hardening  model.  The  reasons  for  introducing  the  strain 
softening  model  in  this  thesis  are  as  follow. 

1.  Part  of  the  strain  softening  model  is  used  in  Chapter  5 
to  improve  the  capabilities  of  the  work  hardening  model. 

2.  The  suggestions  made  for  the  use  of  the  hardening  law 
given  in  Section  3.3  are  also  valid  for  the  strain 
softening  model.  It  is  believed  that  a  critical 
evaluation  of  the  hardening-softening  law  of  the  strain 
softening  model  will  be  useful  for  future  developments 
in  the  subject. 

In  this  chapter  the  strain  softening  model  is  presented 
first  as  developed  by  Lade  (1977).  Then,  by  following  a 
similar  line  of  thought  to  that  expressed  in  Chapter  3  , the 

generality  of  the  hardening-softening  law  is  evaluated. 
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Essential  conditions  for  strain-softening  to  take  place  in 
cohesionless  soils  are  also  examined,  and  reference  will  be 
made  to  this  when  the  finite  element  results  of  passive 
earth  pressure  tests  are  evaluated  in  Chapter  7. 


4.2  STRESS-STRAIN  RELATIONSHIP 

In  the  strain  softening  model,  it  was  assumed  that  the 


total  strain  increments 
elastic 

plastic  expansive  part, 
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a  plastic  collapse,  ^d€ 
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Eqn .  4.1 


The  parts  of  the  total  strain  are  illustrated 
schematically  in  Figure  4.1.  The  calculation  of  each  part  is 
discussed  in  the  following. 


4.2.1  Elastic  Strain  Increments 


d€* 


,  are  calculated 


The  elastic  strain  increments, 
using  Hooke's  law.  The  unloading-reloading  branch  of  a 
stress-strain  curve  is  used  to  determine  the  elastic  modulus 
as  presented  for  the  work  hardening  model  in  section  2.2.1. 
Poisson's  ratio  is  chosen  equal  to  0.2  based  on  the  work  by 
Duncan  and  Chang  (1970),  Calladine  (1973). 
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4.2.2  Plastic  Strain  Increments 

Plastic  strain  increments  are  divided  into  two  parts: 
plastic  collapse  strains  and  plastic  expansive  strains. 
These  parts  are  calculated  from  two  distinct  stress-s trai n 
relationships.  The  yield  criterion,  flow  rule  and  hardening 
law  used  for  each  part  are  formulated  separately.  The  shape 
of  the  yield  surfaces  related  to  different  parts  of  plastic 
straining  are  illustrated  in  Figure  4.2. 

4.2.3  Plastic  Col  1  apse  Strains 

Contrary  to  the  assumption  made  for  metals  in  the 
classical  theory  of  plasticity,  soils  behave  in  such  a  way 
that  part  of  the  strain  increment  due  to  hydrostatic 
compression  is  irrecoverable.  Similarly,  under  proportional 
loading  part  of  the  strain  increment  is  plastic  in  nature. 
It  is  believed  that  the  plastic  collapse  strains  are 
produced  by  re-arrangement  of  the  grain  structure  and  this 
results  in  a  volumetric  reduction.  The  development  of  the 
stress-s trai n  relationship  is  as  follows. 


1 .  Yield  Cr i ter  ion 

The  cap  yield  surface  used  by  Lade  (1975)  is  a 
sphere  with  the  center  in  the  origin  of  the  principal 
stress  space  as  shown  in  Figure  4.2,  and  is  described  by 
the  function 
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fcot=  if  +  2xI2  .  Eqn.  4.2 

where  I-,  and  I2  are  the  first  and  second  stress 
invariants  respectively. 

If  a  stress  increment  results  in  an  increase  in  the 
value  of  the  yield  function  fcol  .  i  .  e .  ,  dfco!  >  0 ,  then 
the  soil  will  undergo  e 1 asto-p 1  as t i c  deformation  while 
work  hardening  takes  place.  This  type  of  yielding  does 
not  result  in  eventual  failure.  The  type  of  yielding 
which  causes  the  conical  yield  surface  to  expand  is 
responsible  for  soil  failure. 


2 .  Flow  Rule 

An  associated  flow  rule  is  used  for  the  calculation 
of  plastic  collapse  strain  increments.  Hence,  the 
plastic  potential  function  is  identical  to  the  yield 
function  and  is  expressed  as: 


floor 


2x1 


Eqn .  4.3 


The  relation  between  the  plastic  collapse  strain 
increments  and  the  plastic  potential  function  is  given 


as : 
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Eqn .  4.4 


A  col  is  the  proportionality  constant  and  its  value 
is  determined  by  the  work  hardening  law. 


3.  Work  Hardening  Law 

The  work  hardening  law,  requried  for  the 
calculation  of  the  magnitudes  of  the  plastic  strain 
increments,  is  expressed  as  an  experimentally  determined 
relation  between  the  total  plastic  work,  WCO|  due  to 
collapse  strains,  and  the  value  of  the  yield  function 

fco,  • 

Wcol  =  C  x  p  *  (“31 


Eqn.  4.5 


where  p  is  the  atmospheric  pressure,  and  C  and  p  are 
a 

constants  which  can  be  determined  from  a  (WCO|/p  ) 

versus  (tCO|/n2)  plot  in  log- log  scale,  as  shown  in 

'a 

Figure  4.3.  It  is  assumed  that  the  work  hardening 
relationship  is  independent  of  the  stress  path. 

The  value  of  Acoi  can  be  determined  from  the 


following  expression. 
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where 


Eqn .  4.6 


— . ^  ‘'7 

The  derivation  of  Eqn.  4.6  is  given  by  Lade  (1975) 
and  will  not  be  repeated  here. 

By  substituting  Eqn.  4.6  into  Eqn.  4.4  and 
calculating  the  partial  derivatives  of  the  plastic 
potential  function  ,  the  final  form  of  the  stress-strain 
relationship  for  plastic  collapse  strains  becomes 
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Eqn.  4.8 


where  dWCOi  is  given  by  Eqn.  4.7. 


4.2.4  P 1 ast ic  Expansive  Strains 

This  part  of  the  strain  increment  is  calculated 
according  to  a  stress-strai n  theory  which  resembles  the 
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theory  presented  previously  for  the  work  hardening  law. 

1 .  Yield  Surface 

The  yield  surface  is  assumed  to  have  the  same  shape 
as  the  failure  surface  and  is  expressed  by  Eqn .  4.9. 


Eqn .  4.9 


where  fGXp  is  the  stress  level 

m  is  a  material  contant  which  is  responsible  for 
the  curvature  of  the  yield  surface. 

I-i  ,  I3  are  the  stress  invariants. 

It  should  be  noted  that  the  expression  for  stress 
level  given  here  is  different  from  that  of  Eqn.  2.4  .  In 
the  work  hardening  model  the  yield  surface  does  not  have 
curva ture . 

As  the  value  of  fexp  increases,  the  yield  surface 
expands  symmetrically  around  the  hydrostatic  axis.  When 
the  stress  level  is  increased  to  its  peak  value,  i.e., 
fGxp=  n1  the  failure  state  is  reached.  Beyond  this 
point  soil  can  not  support  any  further  increase  in 
stress  level,  and  is  not  capable  of  sustaining  previous 
stress  levels  under  continued  straining. 

The  parameters  n1  and  m  can  be  found  by  plotting 
-27)  versus  (p  7I-,  )  at  failure  in  a  log-log 
scale.  In  Figure  4.4,  n-,  is  shown  as  the  intercept  of 
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the  straight  line  with  (PQ/I1)=1  and  m  as  its  slope. 


The  effect  of  (J2  is  included  in  this  failure 
cr i ter  ion . 

2.  Flow  Rule 

A  non- assoc i a  ted  flow  rule  is  employed  in  the 
development  of  the  stress-strai n  relationship  to 
calculate  the  plastic  expansive  strain  increments. 
Consequently,  the  plastic  potential  function  which  is 
given  by  Eqn .  4.10  is  different  from  the  yield  function, 


Eqn .  4.10 


where  gQV_  is  the  value  of  the  plastic  potential  and 

V  a  P 

n2  is  a  constant  for  given  values  of  fGXp  and  (J3  . 
The  value  of  n2  can  be  found  by  using  the  following 
expression . 


Eqn.  4.11 


in  which  S,  R  and  t  are  material  contants. 

The  stress-strain  relationship  for  plastic 


expansive  strains  is  then  derived  according  to  Eqn.  4.12 
in  incremental  form. 
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d  £exp 
»J 

The  magnitude  of 
harden i ng- soften i ng 


determi ned 


1  aw . 


Eqn . 


from 


4.  12 


the 


3.  Hardening-Softening  Law 

An  experimentally  determined  relation  between 
plastic  work  and  stress  level  is  used  as  the 
hardening-softening  law  for  the  strain  softening  model. 
The  term  stress  level  is  defined  in  Eqn.  4.9.  The 
plastic  work  is  calculated  from 


W 


exp 


=  J°r 


d€- 


exp 


•J 


Eqn.  4.13 


where  Gy  d  £y 


exp 


is  the  plastic  work  done  per  unit 
volume  during  a  strain  increment  d  £yXp 

Figure  4.5  shows  the  variation  of  the  total  plastic 
work  ,  Wexp  ,  with  the  stress  level,  fexp  ,  and  the 
confining  pressure,  G^  ,  as  obtained  from  triaxial 
tests.  It  is  assumed  that  the  peaks  of  all  fexp  versus 
Wexp  curves  occur  at  fexp  =  n-|  .  The  amount  of  plastic 
work  required  to  reach  the  peaks  is  dependent  on  the 
confining  pressure.  Since  there  are  two  Wexp  values  for 
a  single  fexp  on  each  curve,  the  value  of  Wexp  at  the 
peak  stress  level  is  used  to  distinguish  between  the 
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points  on  the  ascending  and  descending  parts. 

Experimental  results  shown  on  Figure  4.5  are 
approximated  by  exponential  functions  and  the  following 
expression  gives  this  relationship. 

j_ 

f exp=  a  X  e  bWex?  X  .  Eqn .  4.14 


The  parameters  a,  b  and  q  are  constants  for  a  given 
value  of  the  confining  pressure  (J3  .  As  a  limitation 
to  Eqn.  4.14,  Lade  (1975)  indicated  that  the  value  of 
fexp  decreases  asymptotically  to  zero  for  very  large 
values  of  Wexp  • 

The  variation  of  the  total  plastic  work  at  peak 
stress  level  is  approximated  by  Eqn.  4.15  as  a  function 
of  confining  pressure  (J3  . 


w 


peak 


Eqn .  4.15 


where  P  and  1  are  constants  and 

p  is  the  atmospheric  pressure, 
a 

The  values  of  q,  a  and  b  can  be  calculated  from  the 
following  expressions: 


q  =  GU  +  (3 


Eqn.  4.16 
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where  Qt  and  /3  are  constants, 


a  =  n 


exp 


Eqn .  4.17 


w 


peak 


where  e  is  the  base  for  natural  logarithms  and 
n1  is  the  value  of  peak  stress  level. 


b  = 


1 


qxW 


Eqn.  4.18 


peak 


Based  on  Eqn.  4.14,  the  increment  of  plastic  work 
is  expressed  as 


Eqn.  4.19 


The  expression 

for 

the 

eva 1 uat i on 

of 

the 

propor t i ona 1 i ty  constant 

^  exp 

is  given  by 

Eqn . 

4.20 

and  its  derivation  can 

be 

found  i n 

Lade  (1975). 

d  WeXp 


m.n2 


Eqn.  4.20 
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4.3  EVALUATION  OF  HARDEN  I NG- SOFTEN  I NG  LAW  FJDR  PLASTIC 
EXPANSIVE  STRAINS 

The  test  results  shown  in  Figure  4.5  indicate  the 
dependence  of  Wexp  versus  fGxp  relationship  on  the 
confining  pressure  G3  .  Accordingly,  Eqn .  4.14  which  is 
the  expression  for  this  relationship  is  valid  only  for  cases 
where  the  confining  pressure  is  Kept  constant. 
Generalization  of  Eqn.  4.14  into  problems  with  varying 
0~3  is  not  possible,  as  can  be  shown  next. 

Referring  to  Figure  4.6,  it  is  assumed  that  the  initial 
state  of  stress  of  a  soil  element  is  on  the  current  yield 
surface  as  indicated  by  point  A,  with  total  plastic  work 
value  equal  to  WA  .A  stress  increment  is  applied  such  that 
the  state  of  stress  moves  to  point  B  without  exceeding  the 
limits  of  the  current  yield  surface.  No  plastic  deformations 
take  place  for  this  stress  increment  and  therefore  the  total 
plastic  work  remains  the  same,  i.e.  WB  =  WA  .  Since  the 

value  of  (J3  has  changed,  the  value  of  the  parameters  a,  b, 
and .q  have  also  changed  due  to  this  stress  increment. 
Looking  at  Figure  4.6b  point  B  i s  on  a  curve  for  which  the 
parameters  are  based  on  the  new  0*3  and  the  stress  level 

fexp  corresponds  to  point  B  of  Figure  4.6a.  The  only  value 
of  the  total  plastic  work  which  can  satisfy  Eqn.  4.14  for 
point  B  is  indicated  by  WBtrue  .  Since  the  plastic  work  value 

for  point  B  attained  at  the  end  of  the  stress  increment 

going  from  A  to  B  as  in  Figure  4.6a  is  not  equal  to 
WBtrue  ,  Eqn.  4.14  is  no  longer  valid.  Therefore,  a  unique 


. 
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relationship  between  fexp  and  Wexp  does  not  exist,  except 
for  tests  with  (J3  constant. 

In  order  to  use  this  hardening-softening  relationship 
established  for  (J3  =  constant  tests  in  the  analysis  of 
general  three  dimensional  problems  an  assumption  has  to  be 
made.  This  assumption  is  as  follows: 


Considering  Figure  4.7,  at  any  stress  level  fexp  ,  a 
tangent  to  the  f0Xp  versus  Wexp  curve  for  CT3  =  constant  at 
a  Known  state  of  stress  provides  the  relationship  between 


dfexp  and 

dwexp 

It 

should 

be 

noted  that 

for  each 

fexp  versus 

Wexp  curve 

there  exist 

two 

points 

with  the 

same  fexp 

value, 

one 

before 

and 

the  other 

after  the 

peak.  Lade's  approach  to  differentiate  between  the  hardening 
and  the  softening  parts  of  these  curves  is  to  compare  the 
current  values  of  Wexp  and  Wpeak  •  If  Wexp  is 
less  than  WpeQk  ,  hardening  is  taking  place  and  therefore 
f ex p  is  still  on  the  ascending  part  of  the  curve.  If 
Wex p  is  greater  than  WpeQk  ,  the  descending  part  of  the 
curve  which  corresponds  to  strain  softening  has  to  be  used. 
This  rule  can  no  longer  be  accepted  valid  within  the  context 
of  the  new  approach  as  explained  next. 

In  general,  Wexp  ,  fexp  and  (J3  are  not  uniquely 
related  to  each  other.  The  total  plastic  work  Wexp  at  any 
point,  for  example  at  point  D  on  Figure  4.6b,  may  have  any 
value  including  the  larger  values  than  WpGak  for  the 
corresponding  confining  pressure,  O3  .  This  situation  may 
arise  for  certain  stress  paths  leading  to  point  D.  As  an 
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example  of  this,  a  stress  path  is  considered  which  results 
in  no  change  to  the  stress  level  feXp  but  a  large  decrease 
in  confining  pressure.  In  Figure  4.6  the  results  of  such  a 
stress  change,  i.e.,  from  point  C  to  D  is  shown.  At  point  C 
the  total  plastic  work  is  WeXp  which  has  a  larger  value 
than  .  During  the  stress  i ncrement , '  the  stress 

point  remains  on  the  current  yield  surface  as  shown  on 
Figure  4.6a.  Therefore,  when  point  D  is  reached,  there  will 
not  be  any  change  in  total  plastic  work  and  W^XD  will  have 

,  c  c  (a3a)  d 

the  same  value  as  Wexp.  Since  wexp >^peak  *  wexp  is  a1so 

(0"Q)  ^  v 

larger  than  WpeQ^  .  It  should  be  noted  that  at  point  D 

the  soi 1  element  has  not  yet  fai led  and  before  the  peak 
stress  level  is  reached  more  hardening  has  to  take  place  and 
more  plastic  work  has  to  be  done.  Therefore,  WpeQk 

values  can  not  give  a  correct  indication  of  whether  the 
ascending  or  descending  section  of  fexp  versus 

WeXp  curves  to  be  used. 

It  seems  that  the  only  way  to  decide  which  part  of  a 
curve  is  to  be  employed  is  to  keep  track  of  the  stress  level 
at  all  times.  Once  the  peak  stress  level  is  attained, 

further  loading  will  always  require  the  strain  softening 
sections  of  fexp  versus  Wexp  curves  to  be  utilized. 
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4.4  REMARKS  ON  HARDENING  LAWS  FOR  SOILS 

One  of  the  hardening  hypotheses  in  plasticity  theory  is 
stated  as:  the  amount  of  hardening  depends  on  the  total 
plastic  work  (Mendelson,  1968).  The  yield  criterion  is 
wr i t  ten  as 


Eqn .  4.21 


with  the  implication  that  the  resistance  to  further  yielding 
depends  only  on  the  amount  of  work  which  has  been  done  on 
the  material.  The  work  hardening  and  also  the  strain 
softening  models  presented  here  adopt  this  hypothesis.  In 
both  models  the  stress  level,  confining  pressure  and  total 
plastic  work  form  a  unique  relationship,  hence  the  plastic 
strain  increments  can  be  found  uniquely  for  a  given  stress 
increment  if  the  state  of  stress  in  soil  is  known.  In  this 
Chapter  and  also  in  Chapter  3  ,  it  has  been  shown  that  such 
a  unique  relationship  for  f,  O3  and  Wp  does  not  exist. 
For  the  same  state  of  stress,  the  total  plastic  work,  Wp  , 
may  assume  almost  any  value.  It  appears  that  the  work 
hardening  hypothesis  used  in  the  classical  theory  of 
plasticity  does  not  apply  to  soils  and  the  subject  requires 
further  study. 

To  secure  the  uniqueness  of  plastic  strain  increments 
an  i nterpretat ion  of  the  hardening  laws  of  Lade  is  proposed 
in  this  thesis.  Although  satisfactory  results  are  obtained, 
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there  is  no  real  reason  why  a  relationship  between  d f  and 
d  Wp  obtained  from  a  conventional  triaxial  test  data 
should  apply  to  strongly  path  dependent  soil  behavior  in  any 
other  loading  conditions.  Yet  it  is  a  common  practice  in 
metal  plasticity  to  use  a  hardening  relationship  obtained 
from  one  type  of  test  to  predict  material  behavior  under 
other  test  conditions.  This  is  in  fact  one  of  the  basic 
assumptions  of  the  classical  plasticity  theory  (see  Hill, 
1950).  Similarly,  by  using  the  same  approach,  good 
predictions  for  laboratory  tests  on  soils  were  obtained  as 
shown  in  Chapter  3.  With  the  interpretations  presented  here 
for  the  hardening  and  softening  laws,  Lade's  models  can  now 
be  employed  for  the  solution  of  boundary  value  problems. 
Consequently,  the  evaluation  of  these  models  can  be 
completed  with  their  application  to  practical  problems. 


4.5  CONDITIONS  ESSENTIAL  FOR  STRAIN  SOFTENING 

Strain  softening  behavior  of  material  may  play  an 
important  role  in  both  movements  and  instability  of  soil 
structures.  The  progressive  loss  of  strength  with  increased 
deformation  may  cause  a  significant  redistribution  of  the 
stress  and  eventually  a  progressive  failure  of  the 
structure.  The  Key  points  which  play  an  important  role  in 
the  development  of  strain  softening  are  discussed  next  by 
using  conventional  triaxial  test  results.  Figure  4.8  and  4.9 
show  s t ress - s t r a i n- vo 1 ume  change  data  for  loose  and  dense 
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sands.  As  illustrated  on  these  figures: 

1.  The  behavior  of  soil  during  the  shear  stage  is  dependent 
on  the  relative  density  of  the  specimen  attained  at  the 
end  of  all-round  pressure  application.  Under  continued 
straining  after  the  peak  resistance  is  reached,  large 
decrease  in  resistance  is  mostly  related  to  dense  sand. 

2.  An  increasing  confining  pressure  decreases  the  tendency 
to  dilate,  increases  the  strain  to  failure  and  reduces 
the  brittle  characteristics  of  the  stress-strain  curves. 
For  high  confining  pressures  even  dense  sand  does  not 
exhibit  any  significant  softening,  and  the  volume 
changes  are  entirely  compressive.  When  sheared  at  low 


pressures 

,  dense  sands  are 

a  1 1 owed 

to 

di late, 

and 

softening 

becomes  apparent. 

On 

the 

other 

hand , 

sands 

wi th  any 

initial  densi ty  can 

be 

sheared 

wi thout 

any 

volume  change  at  failure  by  applying  an  appropriate  cell 
pressure.  This  pressure,  correspondi ng  to  zero  volume 
change,  is  called  the  critical  confining  pressure.  An 
example  of  it  is  given  in  Figure  4.10  for  samples  of 
Sacramento  River  sand  tested  at  different  void  ratios 
(Lee  and  Seed,  1967).  Points  between  the  origin  of  the 
coordinate  axes  and  the  curve  shown  on  this  figure 
indicate  the  state  of  soil  which  is  most  likely  to 
display  strain  softening  characteristics.  It  appears 
then,  only  the  combined  effects  of  relative  density  and 
confining  pressure  can  reveal  the  nature  of  post  peak 


behavior . 
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Evidently,  for  strain  softening  to  take  place,  the 
stress  level  should  go  high  enough  to  reach  the  peak  value. 
Also,  strains  should  be  large  enough  to  pass  well  beyond  the 
values  cor respondi ng  to  the  peak  strength  before  the  effect 
of  softening  becomes  appreciable. 

These  conditions  are  not  satisfied  in  all  boundary 
value  problems.  When  they  do,  they  might  be  confined  to  a 
small  portion  of  the  structure  and  therefore  might  not 
affect  the  overall  behavior.  When  widespread  areas  of  a  soil 
mass  are  affected,  or  the  deformations  are  concentrated  in 
discrete  zones  then  the  problem  requires  a  realistic 
analysis.  The  effects  of  strain  softening  can  be  studied 
only  by  analyzing  the  complete  soil  structure  under  load. 
The  soil  type  and  triaxial  test  data  alone  are  not  enough 
for  making  judgment. 
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Figure  4.1  Schematic 
collapse  and  plastic 
compression  test.  (After 


illustration  of  elastic,  plastic 
expansive  strains  in  triaxial 
Lade  , 1975) 
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Figure  4.2  Location  of  yield  cap  relative  to  conical  yield 
surface  shown  in  triaxial  plane.  (After  Lade  ,1975) 
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Figure  4.3  Relation 
and  the  value  of  fco| 
(After  Lade  ,  1975 ) 


between  plastic  collapse  work,  Wco[ 
for  loose  Sacramento  River  sand 
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Figure  4.4  Determination  of  the  value  of  n^  and  m  involved 
in  the  failure  criterion  for  loose  Sacramento  River  sand. 


( After  Lade  , 1975 ) 


, 


' 


. 


exp 


73 


74 


Figure  4.6  Results  of  a  stress  change  a)  stress  paths 
b)We^p  versus  fexp  relationship. 
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Figure  4.7  Assumption  used  for  d feXp  versus  dWexp 
re  1  at ionshi p . 
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Figure  4.8  Measured  stress-strain  and  volume  change  behavior 
for  loose  Sacramento  River  sand  in  drained  triaxial 

test.  (After  Lee  and  Seed,  1967). 
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Figure  4.9  Measured  stress-strain  and  volume  change  behavior 
for  dense  Sacramento  River  sand  in  drained  triaxial 

test.  (After  Lee  and  Seed,  1967). 
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0  5  10  15  20  25  30 

Critical  Confining  Pressure  ( Kg/cm2  ) 


Figure  4.10  Relationships  between  void  ratio  after 
consolidation  and  critial  confining  pressure.  (After  Lee  and 
Seed ,  1967). 
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CHAPTER  5 


WORK  HARDENING  MODEL  WITH  A  CAP 

5. 1  INTRODUCTION 

Lade's  work  hardening  model  has  a  number  of 
limitations.  As  discussed  in  Section  2.3,  one  of  these 
limitations  is  unsatisfactory  modelling  in  proportional 
loading.  Similarly,  for  all  stress  paths  shown  in  Figure 
5.1,  the  model  predicts  only  recoverable  strains  which  is  in 
contradiction  with  observed  soil  behavior.  The  cone  shaped 
yield  surface  of  the  model  without  a  cap  is  shown  in  Figure 
5.2a.  Any  stress  increment  along  a  stress  path  remaining 
within  that  cone  does  not  produce  any  plastic  strains.  In 
reality,  plastic  strains  may  take  place  for  some  stress 
paths  in  that  zone.  It  is  proposed  here  to  overcome  this 
limitation  alone  by  attaching  a  cap  type  yield  surface  to 
the  open  end  of  the  conical  yield  surface.  The  details  of 
the  type  of  cap  used  and  the  improvements  obtained  in 
predictions  are  presented  in  this  Chapter. 

5.2  TYPE  OF  CAP 

In  order  to  model  soil  behavior  more  correctly  cap  type 
yield  surfaces  have  been  proposed  in  the  past  by  many 
researchers  starting  with  Drucker  et  al . ,  (1957).  Depending 
on  the  model,  a  spherical  or  some  other  convex  shape  yield 
surface  is  employed  as  a  cap.  If  the  state  of  stress  changes 


79 


' 

. 

, 


I  • 


. :  / 11  *rf  <fG*i  MfHPf 


80 


such  that  this  cap  is  pushed  out  as  in  Figure  5.2b,  then 
plastic  strains  are  predicted.  Among  several  cap  type  yield 
surfaces  available  in  the  literature,  the  cap  yield  surface 
used  by  Lade  (1975)  as  part  of  his  strain  softening  model 
was  selected  as  a  cap  for  the  present  work,  and  is  used  in 
connection  with  his  work  hardening  model.  The  corresponding 
flow  rule  and  the  work  hardening  law  are  also  taken  from  his 
work  without  any  change.  The  details  of  the  yield  criterion 
are  given  in  Section  4.2.3. 

5.3  UTILIZATION  OF  CAP  TYPE  YIELD  CRITERION 

Although  the  cap  yield  surface  of  Lade's  strain 
softening  model  was  chosen  without  any  change  in  its 
formulation,  as  the  desired  cap  for  the  work  hardening  model 
,  the  conditions  under  which  it  is  used  are  not  the  same  in 
the  present  development.  In  order  to  clarify  the  differences 
in  using  the  same  cap  yield  criterion  for  these  two  models, 
it  is  necessary  first  to  elaborate  on  the  corresponding 
sections  of  the  strain  softening  model. 

In  the  strain  softening  model  total  strain  increments 
are  calculated  as  the  sum  of  elastic,  plastic  collapse  and 
plastic  expansive  parts.  It  is  assumed  that  each  part  can  be 
calculated  separately  and  then  superposed.  An  important 
point  here  is  that  any  stress  change  causing  the  yield  cap 


to  expand,  see  Figure 


without  any  condition  on  where  the  conical  yield  surface  is 
located  at  that  moment.  For  example  the  stress  paths  from  A 
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to  B  and  from  P  to  D  of  Figure  5.3a  and  5.3b  respectively, 
require  that  the  plastic  collapse  strains  be  calculated  and 
superposed  on  the  other  parts  of  strain  as  indicated  by  Eqn. 
4.1.  The  significance  of  the  stress  path  from  A  to  B  i s  that 
it  causes  only  the  cap  yield  surface  to  expand  while  the 
conical  yield  surface  stays  where  it  was  before.  On  the 
other  hand  the  stress  path  from  P  to  D  results  in  expansions 
of  both  yield  surfaces. 

Eqn.  2.1  which  is  somewhat  similar  to  Eqn.  4.1 

represents  the  concept  of  dividing  the  total  strain 

increments  into  elastic  and  plastic  parts  for  the  work 

hardening  model.  However,  the  plastic  strain  increments, 
Pi 


{d£: 


,  are  not  further  sub-divided  into  plastic  collapse 
and  plastic  expansive  parts  as  in  the  case  of  the  strain 
softening  model.  For  example,  the  plastic  strain  increment 
taking  place  due  to  the  stress  change  from  P  to  D  in  Figure 


5.2c,  is  represented  as 


HI 


alone  by  using  the  work 


hardening  model  while  the  same  plastic  strain  increment  is 
calculated  as  the  sum  of  two  components,  i.e., 


exp 

'j 


when  using  the  strain  softening  model.  For 


proportional  loading,  such  as  the  stress  path  1  of  Figure 


5.1,  \  d£y 


becomes  zero  while  both 


<1 


and 


<ixp 


have  non-zero  values.  It  is  clear  thatjd6jj|  of  Eqn.  2.1 

f  €XP]  ^ 

does  not  have  the  same  meaning  as  sd£jj  j  of  Eqn.  4.1. 

When  the  stress  changes  take  place  along  any  one  of  the 
stress  paths  shown  in  Figure  5.1,  the  work  hardening  model 
without  a  cap  predicts  only  elastic  strains.  On  the  other 
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hand,  the  work  hardening  model  with  a  cap  is  able  to  predict 
plastic  collapse  strains  in  addition  to  the  elastic 
increments  of  strain  for  the  same  stress  paths. 

A  summary  of  requirements  for  calculating  different 
parts  of  total  strain  increments,  by  using  the  work 
hardening  model  with  a  cap,  is  now  provided. 

a)When  a  stress  path  remains  within  the  current  conical 
yield  surface  (indicated  as  the  shaded  area  in  Figure  5.2a) 
the  total  strain  increments  are  calculated  by 


* 
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Eqn .  5.1 


Calculation  of  |d€jj°  is  given  by  Eqn.  4.8,  and  the 
proportional  loading  path  is  included  in  its  zone  of 
application.  Of  course,  df  has  to  be  bigger  than  zero  for 


d€ 


col 


to  be  calculated,  otherwise  it  must  be  zero. 


b ) I f  the  stress  increment  is  along  a  stress  path  directed 
outward  from  the  conical  yield  surface,  as  in  Figure  5.2c, 
and  5. 2d,  the  total  strain  increments  are  calculated  by 


d€ij 


Hd4N 


Eqn .  5.2 


which  is  identical  to  Eqn.  2.1. 

c ) When  the  state  of  stress  follows  the  path  from  P  to  D  as 
shown  in  Figure  5.2c,  the  plastic  collapse  strains  are  not 


calculated  separately  but  • 


d€jU  of  Eqn.  2.1  represents  the 
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sum  of  all  parts  of  plastic  strain  increments.  Care  has  to 


be  taken  for  this  type  of  stress  path.  Although 


1  s 


not  calculated,  the  value  of  fcol  has  to  be  updated.  In  other 
words,  the  cap  yield  surface  has  to  be  moved  to  a  new 
location  passing  through  point  D. 

5.4  JUSTIFICATION  FOR  WORK  HARDENING  MODEL  WITH  A  CAP 

As  presented  earlier  in  section  4.1,  a  cap  improves  the 
capabilities  of  the  work  hardening  model  and  the  calculation 
of  plastic  strains  which  may  take  place  due  to  stress 
increments  along  certain  paths  remaining  within  the  conical 
yield  surface  becomes  possible.  This  point  has  been 
emphasized  enough  already. 

Comparisons  of  observed  stress  strain  curves  with  those 
predicted  by  the  work  hardening  model  with  a  cap  are 
presented  in  Section  5.5.  to  illustrate  the  possible 
improvements  that  can  be  obtained  in  prediction  of 
laboratory  test  results. 

The  significance  of  a  cap  yield  surface  in  the  analysis 
of  a  geotechnical  engineering  problem  may  be  predicted, 
provided  that  the  stress  paths  followed  by  the  soil  elements 
comprising  the  structure  are  known.  If  the  majority  of  the 
elements  undergo  proportional  loading  while  the  remaining 
elements  experience  a  very  small  amount  of  plastic  expansive 
straining,  the  yield  cap  will  have  a  dominant  role  in  the 
prediction  of  behavior.  Without  a  cap,  the  work  hardening 
model  predicts  mostly  elastic  strains  for  such  problems. 
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There  are  important  soil  structures  where  stress  paths 
generally  follow  proportional  loading,  for  example  earth 
dams  as  indicated  by  Eisenstein  and  Law  (1979). 

For  the  justification  to  be  complete,  the  following 
question  is  to  be  answered:  why  is  there  a  need  for  a  work 
hardening  model  with  a  cap  while  there  already  exists  a  more 
sophisticated  model  which  has  a  cap  and  can  take  the  strain 
softening  into  account?  There  are  a  number  of  reasons  for 
this  effort  as  discussed  next. 

1.  Not  all  soils  under  all  circumstances  exhibit  strain 
softening.  Key  conditions  for  strain  softening  to  take 
place  are  presented  in  Chapter  4.  Therefore  there  may 
not  be  a  need  for  more  sophistication  and  the  work 
hardening  model  with  a  cap  may  be  just  as  good. 

2.  For  a  problem  where  strain  softening  does  not  occur,  an 
unnecessarily  refined  analysis  using  the  strain 
softening  model  will  be  too  costly. 

3.  The  calculations  required  to  obtain  the  parameters  for 
the  softening  model  are  more  involved. 

4.  The  development  of  a  finite  element  program  for  a  strain 
softening  model  is  more  difficult. 

5.  One  way  of  determining  the  significance  of  modelling  the 
strain  softening  behavior  of  soil  in  analysis  of  soil 
structures  ,  is  to  compare  the  predictions  of  two  finite 
element  programs  developed  for  constitutive  laws  with 
different  capabilities  for  modelling  strain  softening 
alone  but  otherwise  exactly  the  same.  It  is  believed 
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that  the  Lade' s  strain  softening  model  and  the  work 
hardening  model  with  a  cap  can  be  used  for  such  a 
comp a r i son . 

6.  Finally,  the  importance  of  using  a  cap  yield  criterion 
in  the  analysis  of  boundary  value  problems  can  be 
investigated  by  examining  the  differences  of  predictions 
made  by  two  finite  element  programs  developed  for  models 
such  as  Lade's  work  hardening  model  with  and  without  a 
cap  yield  surface. 


5.5  PREDICTIONS  OF  STRESS-STRAIN  BEHAVIOR 


The  measured  behavior  of  two  different  soils  tested 
along  various  stress  paths  are  presented  next  with  the 
predictions  of  the  work  hardening  model  with  a  cap. 

The  first  set  of  test  results,  using  Ottawa  sand  was 
provided  by  the  organizing  committee  of  the  NSF/NSERC  North 
American  Workshop  on  "Plasticity  Theories  and  Generalized 
Stress- Stra i n  Modelling  of  Soils"  held  in  McGill  University, 
Montreal  in  May  1980.  The  parameters  required  in  the  present 
work  for  the  calculation  of  the  elastic  and  plastic  collapse 
strains  were  obtained  from  conventional  triaxial  compression 
and  hydrostatic  compression  test  results.  The  stress-strai n 
curve  as  predicted  and  the  experimental  results  for  the 
stress  path  along  the  hydrostatic  stress  axis  are  plotted  in 
Figure  5.4.  The  results  show  that  there  is  a  significant 
improvement  in  predictions. 

The  second  experimental  data  is  taken  from  the 
published  work  of  Lade  and  Duncan  (1976),  in  which  the 
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results  of  tests  on  loose  Monterey  No.O  sand  were  provided 
for  a  stress  path  representing  proportional  loading.  Those 
measured  results  and  the  preditions  made  using  the  cap  model 
are  presented  in  Figure  5.5.  The  value  of  Poisson's  ratio 
used  in  this  calculation  is  0.26.  Again,  a  comparison 
suggests  that  by  adding  a  cap  yield  criterion,  some 
improvement  can  be  achieved  over  the  work  hardening  model 
wi thout  a  cap . 

The  results  of  experiments  on  loose  Monterey  No.O  sand 
tested  along  a  stress  path  with  increasing  deviator  stress 
and  constant  confining  pressure  after  previous  loading  to  a 
higher  stress  level  were  also  given  by  Lade  and  Duncan 
(1976).  These  test  results  and  the  predictions  made  using 
the  present  model  for  the  soil  behavior  along  this  stress 
path,  are  shown  on  Figure  5.6.  Although  a  slight  improvement 
is  achieved,  especially  in  the  predictions  of  volumetric 
strains,  there  seems  to  be  more  room  for  further  improvement 
in  modelling  soil  behavior  for  this  kind  of  stress  path  and 
stress  history.  It  should  be  noted  that  the  stress  path  from 
point  A  to  point  B  involves  reloading  because  the  soil 
element  was  subjected  to  an  earlier  higher  stress  level 
before  point  A  was  reached.  Therefore  only  elastic  strains 
are  predicted  for  the  section  between  point  A  and  point  P 
when  the  loading  follows  the  path  from  A  to  B.  The 
s t ress - s t r a i n  curves  in  Figure  5.6  show  that  the  largest 
difference  between  predicted  and  measured  curves  is  at  point 
P.  This  difference  does  not  change  to  a  noticeable  degree 
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under  subsequent  loading.  It  appears  that  during  large 
unloading-reloading  cycles  soils  tend  to  "forget"  the 
maximum  stress  level  attained  previously.  This  contradicts 
an  assumption  used  in  the  development  of  Lade's  work 
hardening  and  strain  softening  models.  In  reality,  the 
current  yield  surface  during  unloading-reloading  cycles  may 
not  stay  at  the  level  corresponding  to  the  maximum  stress 
level  reached  in  its  history,  but  may  shrink  or  shift  its 
position  or  even  may  take  another  shape. 


5.6  CONCLUSIONS 

In  this  Chapter  a  cap  type  yield  criterion  was 
incorporated  into  Lade's  work  hardening  model.  The 
improvement  that  can  be  obtained  by  using  a  cap  yield 
criterion  was  demonstrated  by  comparing  experimental  results 
with  predictions  made  by  the  models  with  and  without  such  a 
cap.  In  all  cases  examined,  the  predictions  made  by  adding  a 
cap  to  the  model  showed  a  better  agreement  with  the 
experimental  data.  The  degree  of  success  attained  in 
modelling  soil  behavior  was  dependent  on  the  stress  path  and 
stress  history.  For  certain  stress  paths  involving 
unloading-reloading,  improvement  was  not  significant.  Cyclic 
loading  requires  more  attention. 

The  work  hardening  model  with  a  cap  links  the  work 
hardening  model  without  a  cap  and  the  strain  softening 
model,  and  altogether  they  provide  a  selection  of  different 
models  with  improved  capabilities. 
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Furthermore,  the  procedure  proposed  in  Section  5.2  to 
combine  two  models  may  be  used  by  other  model  developers. 
Plastic  strain  increments  need  not  be  divided  into  plastic 
collapse  and  plastic  expansive  parts  for  all  stress  paths. 
If  such  a  division  exists,  then  certain  assumptions  have  to 
be  made  in  dividing  the  total  plastic  strain  increment  into 
its  parts.  For  a  general  stress  increment,  direct 
measurement  of  separate  plastic  collapse  and  plastic 
expansive  strain  increments  is  not  possible. 

Finally,  the  finite  element  procedures  become  easier 
for  constitutive  laws  dealing  with  smaller  number  of 
stress-strai n  relationships  within 


the  model . 
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Figure  5 .  1 
i ncrements 


Example  of  stress  paths  where  no  plastic  strain 
are  predicted  by  the  work  hardening  model  without 


a  cap. 
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Figure  5.2  The  effect  of  different  stress  paths  on  the 
expansion  of  yield  surfaces  for  work  hardening  model  (a) 
without  a  cap  (b,  c,  d)  with  a  cap. 
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Figure  5.4  Hydrostatic  compression  of  Ottawa  sand. 
Predictions  of  work  hardening  model  with  and  without  a  cap 
yield  criterion.  (  Experimental  data  from  N5F/NSERC  Workshop 
on  Soil  Plasticity,  Montreal,  May  1980) 
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Figure  5.5  Proportional  loading.  Predictions  of  work 
hardening  model  with  and  without  a  cap.  Loose  Monterey  No.  0 
sand.  Experimental  data  after  Lade  and  Duncan  (1976) 
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Figure  5.6  Stress  path  with  increasing  deviator  stress  and 
constant  C3  after  previous  loading  to  a  higher  stress 
level  (Experimental  data  after  Lade  and  Duncan,  1976) 


CHAPTER  6 


FINITE  ELEMENT  PROCEDURES  FOR  WORK  HARDENING  MODEL 
6. 1  INTRODUCTION 

In  order  to  analyze  plane  strain  problems,  a  finite 
element  program  was  developed  by  using  the  displacement 
formulation.  The  main  intention  was  to  examine  the 
capabilities  of  Lade's  work  hardening  model  (without  a  cap) 
to  predict  soil  behavior  under  complex  boundary  conditions. 
Also,  the  difficulties  related  to  the  use  of  the  model  in  a 
finite  element  analysis  were  investigated.  Efficiency, 
simulation  of  construction  stages  or  modeling  soil  structure 
interface  behavior  were  not  the  prime  concern;  their  study 
was  left  until  after  the  usefulness  of  the  model  was  proven. 

The  work  hardening  law  of  the  model  was  used  in  its 
re-evaluated  form  as  presented  in  Chapter  3.  This  led  to  a 
different  constitutive  matrix  for  the  finite  element 
formulation  than  the  ones  proposed  by  Ozawa  and  Duncan 
(1976a),  Duncan  et  al  (1977)  and  Wong  (1978).  Therefore,  the 
results  of  the  finite  element  analysis  of  the  present  work 
and  the  conclusions  reached  are  different  from  the 
calculations  and  arguments  of  all  the  other  studies,  as  will 
be  discussed  in  the  next  two  Chapters. 

A  linear  elastic  finite  element  program  coded  by  Murray 
(1974)  was  used  as  a  starting  point.  Major  changes  were  made 
in  order  to  assign  initial  stresses  to  elements,  to  solve  a 
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system  of  algebraic  linear  simultaneous  equations  having  an 
unsymmetric  coefficient  matrix,  and  to  deal  with  material 
nonlinearity.  Isoparametric  elements  were  employed  in  this 
program. 


6.2  SOLUTION  METHOD  FOR  MATERIAL  NONLINEARITY 

Convergence  problems  encountered  by  others,  such  as 
Naylor  (1975)  and  Wong  (1978)  were  taken  as  a  lesson  in  the 
present  work  for  the  selection  of  a  solution  method  to 
handle  the  material  nonlinearity.  The  solution  technique 
chosen  is  called  the  tangent  stiffness  method.  The  basis  of 
this  procedure  is  the  application  of  the  load  in  small 
increments.  A  new  stiffness  matrix,  reflecting  the  changing 
stress-s trai n  relationship,  is  formed  at  each  increment. 
Thus,  the  procedure  approximates  the  nonlinear  problem  as  a 
series  of  linear  problems. 

To  improve  the  accuracy  of  the  incremental  procedure 
and  to  reduce  the  deviation  from  the  true  solution,  an 
iterative  technique  was  introduced  for  the  tangent  stiffness 
method.  In  Figure  6.1,  the  technique  is  illustrated  for  a 
single  degree  of  freedom  system.  The  starting  point  for  the 
load  increment,  A  R  ,  is  represented  by  point  A  at  which  the 
stress,  CTa  ,  is  known.  The  tangent  stiffness  based  on  the 
state  of  stress  at  point  A  ,  K1=K(0"A  ),  is  calculated  first. 
The  solution  using  the  stiffness  K1  yields  displacement 
i ncrement  AS..  From  this  a  strain  increment  A€1  is 
obtained.  The  stress  increment  Ad,  correspondi ng  to 
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is  then  evaluated  by  using  the  stress-strain  relationship 
formed  for  the  stress  (JA 

The  procedure  described  so  far  constitutes  the  simple 
incremental  tangent  stiffness  method.  The  iteration  scheme 
introduced  here  to  improve  the  approximation  for  the  same 
load  increment  is  as  follows: 

After  the  stress  increment  AO}  is  calculated  (as 
described  above)  a  new  stiffness  K2  is  formed  as  a 
function  of  the  stress  C^+y(  A(Jj  )  ,  and  a  new  estimate 
for  the  displacement  increment  A§2  is  calculated  by 
solving  the  system  of  equations  for  the  same  load 
increment.  Then  the  corresponding  strain  increment  A62 
is  obtained  from  a62  .  Calculation  of  AO}?  is  based 
on  the  stress-strain  relationship  produced  for  the 
stress  C^+y-(  ACT  )  . 

The  iteration  process  can  be  repeated  until  the 
difference  between  two  successive  solutions  is 
sufficiently  small.  In  the  present  work  the  iterations 
were  carried  out  once  or  twice. 


6.3  DERIVATION  OF  CONSTITUTIVE  MATRIX 

The  nonlinear  material  behavior  is  approximated  by 
changes  in  stress-strai n  relationship  during  successive  load 
increments.  Therefore  the  constitutive  matrix  suitable  for 
incremental  finite  element  analysis  has  to  be  formulated 
such  that  the  increments  of  stress  can  be  calculated  for 
given  strain  increments  at  known  stress  state  of  soil. 
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The  following  derivation  is  similar  to  the  one  given  by 
Zienkiewicz  (1971),  with  the  exception  that  the  necessary 
changes  were  introduced  to  take  into  account  the 
non- associ ated  flow  rule  used  in  the  present  study. 

During  an  infinitesimal  increment  of  stress,  the  strain 
increments  taking  place  are  assumed  to  be  separable  into 
elastic  and  plastic  parts. 


d 
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+  d 


Eqn .  6.1 


The  elastic 

strain 

increments  are  related 

to  the 

stress 

i ncrements 

by 

a  symmetric  matrix, 

_ce_ 

,  as 

wr i t  ten 

expl ici t ly 

in  Eqn. 

CO 

CM 

Plastic  strain  increments  were  given 

in  Eqn.  2.7 

.  Thus, 

Eqn . 

6.1  can  be  written  as 

mm 

-1 

l  ,  0g 

J  0  |  O’ 


The  yield  function  may  be  expressed  for  a  work  hardening 
material  as 


Eqn .  6.3 


When  differentiated,  Eqn.  6.3  gives 
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Eqn .  6.4 


Eqn .  6.5 


in  which  a  substitution  is  made  for 


Eqn .  6.6 


Eqn.  6.2  and  Eqn.  6.5  is  now  written  in  matrix  form  as  given 
next . 
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Eqn.  6.7 
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By  inverting  the  square  matrix  of  Eqn.  6.7,  an  expression  is 
obtained  to  determine  the  stress  increments  in  terms  of 
specified  strain  increments. 
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Eqn .  6.8 


where  [C6^]  i s  the  elasto-plastic  constitutive  matrix  and 
is  formed  by  using  Eqn.  6.9.  (see  the  next  page.) 
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The  el asto-pl ast i c  constitutive  matrix 
replaces  the  elastic  matrix  in  incremental  analysis 


when  plastic  as  well  as  elastic  straining  takes  place.  It 


should  be  noted  that 


is  no  longer  a  symmetric 


matrix  because  the  partial  derivatives  of  the  plastic 
potential  and  the  yield  functions  are  not  the  same  in  Lade's 
work  hardening  model. 

6.4  VERIFICATION  OF  FINITE  ELEMENT  PROGRAM 

In  geotechnical  engineering  practice  today,  the  trend 
to  test  the  correctness  of  newly  developed  finite  element 
programs  for  nonlinear  stress-strai n  relationships  is  to 
compare  the  calculated  values  with  laboratory  test  results. 
Here  an  alternative  and  more  accurate  procedure  is  offered. 
Rather  than  using  experimental  results,  predictions  made  by 
the  model  (such  as  the  ones  given  in  Chapter  3)  are  employed 
for  testing  the  results  of  the  finite  element  program.  The 
advantage  of  using  the  predictions  is  to  separate  the 
capabilities  of  the  model  from  the  correctness  and  the 
accuracy  of  the  finite  element  program.  This  approach  can  be 
used  effectively  for  boundary  value  problems  where  the  state 
of  stress  and  strain  are  uniform.  Three  test  cases  chosen 
for  the  verification  of  the  finite  element  program  of  the 
present  study  fpllow. 
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TEST  CASE  I:  ANALYSES  OF  PLANE  STRAIN  TESTS 

Plane  strain  test  results  shown  in  Figure  6.2  and 
Figure  6.3  for  loose  and  dense  Monterey  No.  0  sand  are  used 
to  illustrate  the  verification  procedure.  The  related 
stress-strain  parameters  as  derived  by  Lade  are  listed  in 
Table  6.1.  By  using  the  constitutive  model  in  a  simple 
stress  analysis,  it  is  possible  to  predict  the  plane  strain 
behavior  of  soil  without  resorting  to  finite  element 
calculations.  For  a  given  constant  cell  pressure  and  an 
increment  of  axial  stress  from  a  Known  state  of  stress,  the 
increment  of  intermediate  principal  stress  can  be  calculated 
in  an  iterative  procedure  which  satisfies  the  plane  strain 
condition  as  defined  by  A  £2=  0  •  The  predictions  made  by 
such  an  analysis  are  shown  on  Figure  6.2  and  Figure  6.3.  It 
should  be  noted  that  the  comparisons  of  predictions  with 
experimental  results  give  simply  an  indication  of  the 
capabilities  of  the  model.  For  the  analysis  of  the  same 
problems  using  the  finite  element  method  only  one  element  is 
required  because  the  state  of  stress  and  strain  are  assumed 
to  be  uniform.  Calculated  stress-strain  and  volume  change 
values  by  finite  element  analysis  are  also  plotted  in  Figure 
6.2  and  Figure  6.3.  The  agreement  between  these  results  and 
the  predictions  of  simple  stress  analysis  is  excellent  while 
they  both  deviate  from  the  experimental  curves.  For  example, 
as  shown  in  Figure  6.3,  the  measured  axial  strain  at 
deviator  stress  6.0  Kg/cm  is  1.25  %.  For  the  same  deviator 
stress,  the  simple  stress  analysis  and  the  finite  element 
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method  give  axial  strains  approximately  5.0  %  which  is  four 
times  larger  than  the  experimental  value.  This  comparison 
clearly  indicates  that  finite  element  programs  should  not  be 
tested  against  experimental  results,  but  the  predictions 
made  by  simple  stress  analyses  should  be  used  for  this 
purpose . 

TEST  CASE  II:  Ko-L0ADING  CONDITION 

K0  -  loading  condition  is  satisfied  when  a  soil  element 


is  subjected 

to 

vertical  load  increments  with 

no 

deformat  ions 

a  1 1  owed  i n 

the  horizontal 

di rect ions . 

The 

behavior  of 

soi  1 

under 

these  boundary 

condi t ions 

has 

significant  importance  in  geotechnical  engineering  because 
the  accumulation  of  sediments  is  assumed  to  take  place  under 
similar  conditions.  Analyses  of  the  behavior  of  soil 
deposits  under  subsequent  changes  in  boundary  conditions 
require  the  knowledge  of  in-si tu  stress  state  prior  to  these 
changes  such  as  an  excavation  or  foundation  loading  etc. 
Soil  behavior  under  K0" loading  conditions  has  been  studied 
extensively  in  the  past.  Although  field  measurements  of  the 
stress  state  of  soil  deposits  have  yet  to  be  perfected,  the 
results  of  laboratory  studies  starting  with  the 
comprehensive  investigation  by  Bishop  (1953)  and 
subsequently  by  Brooker  and  Ireland  (1965)  indicated  that 
the  following  expression  proposed  by  Jaky  for  initial 
loading  is  valid  for  all  practical  purposes. 
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Kc  =  1  -  si n  0* 


Eqn .  6.10 


where  K0denotes  the  coefficient  of  earth  pressure  at  rest 

^"horz. 

°~~CJ 

^vert. 

and  (p '  denotes  the  angle  of  shearing  resistance  in  terms  of 
effective  stress. 

As  a  second  test  for  the  correctness  and  accuracy  of 
the  finite  element  program,  soil  behavior  under  K0- loading 
conditions  is  examined.  While  the  vertical  stress  on  a  soil 
element  is  increased,  the  no  horizontal  strain  increment 
condition  can  be  met  if  the  following  equation  is  satisfied. 

d€h  +  d£eh  =  0  .  Eqn.  6.11 

e  p 

where  d€p  and  d€h  refer  to  horizontal  elastic  and 
plastic  strain  increments  respectively.  Again,  without  using 
the  finite  element  program,  the  behavior  of  an  infinitesimal 
soil  element  under  Ko-loading  condition  is  predicted  first. 
Starting  from  a  very  small  hydrostatic  state  of  stress,  the 
vertical  stress  is  increased  in  steps.  For  each  step  the 
horizontal  stress  increment  which  would  satisfy  the 
condition  given  by  Eqn.  6.11  is  searched.  The  predictions 
made  by  this  analysis  for  both  loose  and  dense  sand  are 
shown  in  Figure  6.5  and  Figure  6.6.  For  the  same  boundary 
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conditions,  the  results  of  finite  element  calculations  are 
also  plotted  on  the  same  Figures.  It  appears  that  for  this 
stress  path  too,  the  finite  element  program  gives  accurate 
resu Its. 

There  are  no  experiments  to  support  the  predictions 
shown  in  these  Figures,  but  the  measured  (jj  values  for 
loose  and  dense  sands  are  given  34.8  and  45  degrees 
respectively  by  Wong  (1978).  Using  Eqn .  6.10,  the 
coefficient  of  earth  pressure  at  rest  are  calculated  as 
Ko=0.43  for  loose  sand  and  Ko=0.29  for  dense  sand.  These 
numbers  are  significantly  different  than  the  predicted 
values  which  are  Ko=0.28  for  loose  sand  and  Ko=0.24  for 
dense  sand. 

Predictions  might  be  improved  by  using  the  work 
hardening  model  with  a  cap  since  the  measured  stresses  in  Kq 
-loading  conditions,  such  as  given  by  Andrawes  and  El-Sohby 
(1973)  indicate  a  stress  path  causing  the  cap  yield  surface 
to  expand. 

TEST  CASE  III:  PREDICTIONS  OF  PASSIVE  EARTH  PRESSURE 

BEHAVIOR 

It  was  assumed  that  three  soil  samples  representative 
of  soil  at  different  depths  behind  a  retaining  wall  were 
subjected  to  passive  earth  pressure  tests  individually. 
Figure  6.6  shows  the  hypothetical  retaining  wall,  the 
sampling  locations  and  the  initial  state  of  stress  for  each 
element.  The  soil  type  chosen  for  this  case  is  dense 
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Monterey  No.  0  sand  and  its  stress-strain  parameters  are 
given  in  Table  7.2.  The  behavior  of  these  samples  under 
plane  strain  boundary  conditions  were  predicted  by  a  stress 
analysis  in  which  (Thorz  increased  in  small  steps  while 
(Jvert  was  Kept  constant.  The  predictions  of  this  stress 
analysis  are  shown  in  Figure  6.7.  Calculated  values  by  the 
finite  element  analysis  are  plotted  on  the  same  Figure.  The 
agreement  once  again  is  excellent  between  the  simple  stress 
analysis  and  the  finite  element  calculations.  A  discussion 
on  these  results  in  terms  of  their  engineering  significance 
will  be  given  in  Chapter  7  after  the  experimental  results  on 
passive  earth  pressure  studies  are  reviewed. 


6.5  CONCLUSIONS 

In  this  Chapter  the  finite  element  procedures  employed 
in  the  computer  program  were  presented  and  a  new  approach 
was  proposed  for  testing  finite  element  programs.  As 
demonstrated  by  three  test  cases,  the  finite  element 
analyses  give  the  same  results  as  the  predictions  made  by 
using  the  model  and  simple  stress  analysis.  More  complex 
studies  by  finite  element  analysis  are  expected  to  give  a 
good  indication  of  overall  capabilities  of  the  model  to 
predict  the  behavior  of  soil  structures. 
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DISPLACEMENT  (§) 


Figure  6.1  Single  degree  of  freedom  representation  of 
i ncrement a  1  - i terat i ve  procedure. 
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Figure  6.2  Results  of  plane  strain  analysis  and  experimental 
data  for  loose  Monterey  No  0  sand.  (Experimental  data  after 
Lade,  1972) 
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Figure  6.3  Plane  strain  analysis  and  experimental  data  for 
dense  Monterey  No  0  sand.  (Experimental  data  after  Lade, 
1972) 
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TABLE  6. 1 


LADE'S  STRESS-STRAIN  PARAMETERS 
FOR  MONTEREY  No.  0  SAND 


Parameter 

D  ense 

Sand 

Loose 

Sand 

Kur 

2300 

1  60  0 

n 

0.8 

0.86 

V 

0 

0 

Ki 

103 

58 

ft 

40 

33 

A  1 

0.44 

0.39 

A  2 

15.12 

16.48 

l 

1.  32 

1.17 

M 

0.000255 

0.00068 

rf 

0.957 

0.970 

1  12 
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DIRECTION  OF  WALL  MOVEMENT 
(NO  ROTATION) 


SAMPLE 

DEPTH 

m 

^vert. 

kg/cm2 

^horz. 
kg /cm2 

A 

0.1524 

0.025832 

0.006975 

B 

1.5240 

0.258320 

0.069750 

C 

15.2400 

2.583200 

0.697500 

Figure  6.6  Boundary  conditions  for  the  hypothetical 
retaining  wall  and  the  assumed  initial  state  of  stress  at 
different  depths  before  wall  movement. 


' 

; 

HORIZONTAL  STRESS  (kg/cm2) 


115 


Figure  6.7  Predictions  of  passive  earth  pressure  for  dense 
Monterey  No.  0  sand  and  finite  element  calculations. 


CHAPTER  7 


FINITE  ELEMENT  ANALYSIS  OF  PASSIVE  EARTH  PRESSURE 
7  .  1  INTRODUCTION 

In  recent  years,  several  "non-classical"  stress-strain 
laws  have  been  proposed  to  model  soil  behavior  more 
realistically.  The  predictions  of  some  of  these  models  show 
remarkable  agreement  with  laboratory  test  results.  Their 
capabilities  were  discussed  in  the  NSF/NSERC  North  American 
Workshop  on  "Plasticity  Theories  and  Generalized 
Stress-Strain  Modelling  of  Soils,"  Montreal,  May  1980.  In 
this  workshop  it  also  became  evident  that  the  number  of 
applications  of  "non-classical"  models  in  engineering 
practice  has  been  very  little  so  far.  From  the  practicing 
engineer's  point  of  view,  the  goal  of  developing  a 
constitutive  law  is  not  fully  accomplished  unless  it  can  be 
used  successfully  and  efficiently  in  the  analysis  of  soil 
structures.  Obviously  there  is  a  considerable  amount  of 
effort  to  develop  finite  element  programs  today.  However, 
considering  the  number  of  researchers  working  on  the  subject 
and  a  rather  long  time  lag  between  the  development  of  a 
model  and  its  first  time  engineering  application, 
substantial  amount  of  difficulties  should  be  anticipated  for 
the  newcomers  to  this  field.  Previous  to  the  present  study, 
there  have  been  two  attempts,  one  by  Ozawa  (1974)  and  the 
other  by  Wong  (1978),  to  use  Lade's  work  hardening  model  in 
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a  finite  element  analysis.  Reference  will  be  made  quite 
frequently  to  their  work  in  this  Chapter.  The  work  presented 
in  the  following  aims  essentially  at  finding  an  answer  to 
one  question  which  is  "Does  Lade's  model  work  in  complex 
boundary  value  problems?".  The  answer  to  this  question  has 
not  been  affirmative  so  far.  Strain  softening  behavior  or 
modelling  the  plastic  collapse  strains  is  not  considered  in 
the  finite  element  analysis  of  the  present  study.  Only  the 
work  hardening  model  without  a  cap  yield  surface  is  used  as 
the  stress-strain  relationship  for  soils. 

The  passive  earth  pressure  problem  has  been  chosen  as 
the  field  of  application  to  test  the  capabilities  of  Lade's 
work  hardening  model.  The  reason  for  this  choice  is  the 
availability  of  numerous  experimental  results  of  model  earth 
pressure  tests  in  the  literature.  The  subject  has  been 
studied  extensively  in  the  past  and  with  the  help  of 
advanced  measurement  techniques  for  stresses  and  strains, 
substantial  amount  of  knowledge  has  accumulated  about  the 
behavior  of  earth  retaining  structures.  Thus,  the  scope  of 
this  Chapter  is  to  find  the  experimental  facts  from 
literature  and  compare  them  with  the  results  of  the  finite 
element  analysis.  No  attempt  is  made  here  to  present  a 
comprehensive  review  of  earth  pressure  theories  and  current 
design  methods  for  earth  retaining  structures.  They  can  be 
found  in  textbooks  as  well  as  proceedings  of  several 
international  and  specialty  conferences  on  the  subject. 

To  proceed  systematically,  this  Chapter  is  divided  into 
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several  major  sections.  In  Section  7.2,  a  brief  literature 
review  reflecting  the  highlights  of  the  experimental  studies 
on  passive  earth  pressure  problems  is  presented  and  the 
factors  which  affect  the  behavior  of  earth  retaining 
structures  are  identified.  The  purpose  is  to  evaluate,  at 
least  qualitatively,  the  results  of  finite  element  studies 
presented  here  and  to  investigate  the  level  of  agreement 
between  predicted  and  expected  behavior.  In  order  to  make 
direct  comparisons  between  predictions  and  experimental 
results  in  evaluating  Lade's  work  hardening  constitutive  law 
,  the  model  studies  by  Wong  (1978)  are  used  as  a  case 
history  in  present  work.  The  details  of  this  experimental 
work  are  given  in  Section  7.3.  Unfortunately,  not  many  case 
histories  presented  in  the  literature  also  provide 
sufficient  data  on  soil  properties  to  derive  the  parameters 
required  for  " non-c  1  ass i ca  1 11  stress-strain  laws.  The  matter 
has  to  be  brought  to  the  attention  of  publishers  in 
geotechnical  engineering  field.  Wong's  finite  element 
analyses  of  the  wall  behavior  i.s  summarized  in  Section  7.4. 
In  Section  7.5  the  predictions  made  by  the  finite  element 
analysis  of  present  study  are  compared  with  Wong's 
experimental  and  analytical  results.  This  comparison 
provides  enough  evidence  for  the  usefulness  of  Lade' s  work 
hardening  model  and  its  superiority  to  the  other 
stress-strain  laws  used  by  Wong  (1978)  to  predict  passive 
earth  pressure  test  results.  The  final  Section  is  devoted  to 
the  i nterpretat ion  of  results  found  in  this  Chapter. 


SIM  rt  ^  «Y 

' 


' 


119 


7.2  EXPERIMENTAL  INVESTIGATION  OF  PASSIVE  EARTH  PRESSURE 
PROBLEM  (LITERATURE  SURVEY) 

One  of  the  earliest  experimental  studies  on  passive 
earth  pressure  problems  was  reported  by  Terzaghi  (1934)  who 
observed  the  importance  of  various  modes  of  wall  movement 
and  the  existence  of  differences  between  the  behavior  of 
dense  and  of  loose  backfill.  In  a  series  of  tests  using  dry 
sand  a  rigid  wall  was  rotated  about  its  toe  in  both  an 
active  and  passive  sense.  Based  on  the  measurements  of 
thrust  acting  on  the  wall,  an  average  earth  pressure  was 
obtained  as  a  function  of  wall  rotation.  Figure  7.1  shows 
the  relation  between  wall  movement  and  mobilized  coefficient 
of  earth  pressure  as  provided  by  Terzaghi  (1954).  The 
important  point  to  be  observed  in  this  Figure  is  that  the 
magnitude  of  earth  pressure  is  a  function  of  the  mode  of 
wall  movement,  amount  of  wall  movement  and  soil  properties. 

Based  on  experiments,  Franzius  (1924),  Tschebotar iof f 
and  Johnson  (1953)  indicated  the  importance  of  wall  friction 
(roughness  of  the  surface)  and  the  boundary  conditions  in 
testing  apparatus.  Tschebotar iof f  (1951)  reported  how 
dredging  or  backfilling  may  affect  the  earth  pressure 
distribution  along  a  flexible  bulkhead.  His  work 
demonstrates  the  significance  of  the  construction  procedures 
and  the  flexibility  of  the  retaining  structure  in  pressure 
distribution.  Work  by  Brinch  Hansen  (1953)  and  Narain  et  al . 
(1969)  are  all  part  of  the  effort  to  increase  our 
understanding  of  passive  earth  pressure  problem. 


i'*im  ■  ■  '  WM  '  •  •  ' 

’ 


I  »  ■  o  a :  .  U  :  I,  -  .  I  o 

■ 


120 


In  a  comprehensive  series  of  tests  in  model  scale,  Rowe 
and  PeaKer  (1965)  investigated  in  detail  the  mobilization  of 
normal  and  shear  stresses  acting  on  a  wall  which  can  be 
translated  in  a  preset  direction  with  the  wall  inclined  to 
the  horizontal  at  any  desired  angle.  Measurements  showed 
that  passive  earth  pressure  values  vary  with  the  orientation 
of  the  wall  and  the  direction  it  is  pushed  against  the  soil. 
In  these  tests  the  mobilized  wall  friction  and  the 
coefficient  of  passive  earth  pressure  Kp  values  did  not 
always  reach  their  peak  values  simultaneously.  For  dense 
sand,  there  was  a  significant  drop  in  Kp  after  a  peak  was 
reached.  The  amount  of  movement  required  to  reach  maximum 
Kp  was  in  the  order  of  15-25%  of  wall  height  for  loose  sand 
and  4-10%  for  dense  sand.  It  is  most  likely  that  the 
magnitude  of  these  movements  is  not  acceptable  in  practice. 
Then,  the  analysis  of  passive  earth  pressure  becomes  a 
deformation  problem.  The  authors  have  concluded  that  "a 
correct  solution  of  earth  pressure  calculation  must  await 
the  establishment  of  the  s tress-strai n . . . 1 aws  for  soils 
subjected  to  any  stress  path" 

One  of  the  most  significant  investigation  in  the 
subject  comes  from  Cambridge  University.  James  and  Bransby 
(1970)  reported  observations  of  normal  and  shear  stress 
distribution  along  a  rotating  wall  as  well  as  the  strains 
within  the  sand  mass.  Shear  strain  contours,  which  were 
found  by  using  an  X-ray  technique,  for  different  degrees  of 
wall  rotation  provided  valuable  information  about  the  actual 
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behavior  of  sand  mass.  Based  on  this  experimental  fact,  they 
then  evaluated  the  validity  of  the  assumptions  made  in 
current  methods  of  earth  pressure  calculations.  Some  of  the 
shear  strain  contours  of  their  work  are  reproduced  here  in 
Figure  7.2  for  an  easy  reference  to  be  used  later  in  this 
Chapter.  Figure  7.3  shows  the  rupture  planes  as  observed  by 
James  and  Bransby  (1970)  on  radiographs  exposed  at  three 
stages  of  a  passive  pressure  test  in  dense  sand.  It  is  clear 
that  the  failed  zone  of  material  becomes  larger  with  the 
increase  in  wall  rotation  in  a  fashion  shown  in  this  Figure. 
In  his  Rankine  lecture,  Roscoe  (1970)  summarized  the  results 
of  passive  earth  pressure  tests  on  model  walls  tested  at  the 
Cambridge  University.  Comparisons  of  the  wall  behavior  for 
different  modes  of  wall  movement  indicate  that  when  the  wall 
rotates  about  its  toe  the  required  angle  of  rotation  for 
peak  resistance  is  about  8  degrees.  If  the  same  wall  is 
rotated  about  the  top,  the  peak  is  attained  at  only  1.4 
degrees  of  wall  rotation.  For  pure  horizontal  translation 
the  displacement  needed  for  maximum  resistance  is  equivalent 
to  2  degrees  of  rotation. 

There  have  been  many  more  experimental  studies  carried 
out  in  the  past.  The  details  of  all  these  are  not  included 
in  this  review.  Nevertheless,  it  is  believed  that  the 
essential  points  of  passive  earth  pressure  problem  in  ideal 
conditions  are  established  sufficiently  for  our  purposes. 
What  is  missing  in  the  experimental  work  presented  here  is 
the  measurement  of  intermediate  principal  stresses. 
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For  large  scale  problems  one  has  to  be  careful  about 
size  effects.  Extrapolation  of  findings  of  small  scale  model 
tests  may  not  be  applicable  directly  to  large  scale 
problems.  There  are  other  factors  such  as  in-si tu  state  of 
stress,  construction  sequences,  temperature  changes,  creep 
effect,  seepage  conditions  and  stresses  caused  by  compaction 
which  may  affect  drastically  the  behavior  of  real  soil 
structures . 

The  most  useful  experimental  work  to  the  present  study 
was  carried  out  by  Wong  (1978),  and  forms  the  basis  of 
comparison  for  the  results  of  a  finite  element  analysis 
which  is  presented  in  Section  7.5.  For  this  purpose  the 
related  part  of  Wong's  work  is  presented  next. 

7.3  EXPERIMENTAL  WORK  BY  WONG 

In  a  series  of  passive  pressure  tests  carried  out  by 
Wong  (1978),  an  instrumented  wall  was  rotated  about  its  toe 
into  a  bed  of  sand  as  illustrated  in  Figure  7.4. 
Measurements  were  taken  for  normal  and  shear  forces  on  four 
panels  of  the  wall  at  various  stages  of  rotation.  The 
outcrops  of  the  slip  planes  on  the  surface  were  mapped  at 
the  end  of  each  test.  Tests  were  conducted  by  using  dense 
(Dr  =  90%)  and  medium  dense  ( D r  =  6 5% )  Monterey  No.O  sand  which 
is  a  uniform,  medium  sand  composed  mainly  of  quartz  and 
feldspar.  Sand  particles  are  subangular  to  subrounded. 
Gradation  curves  for  the  Monterey  No.O  sand  used  in  passive 
pressure  tests  by  Wong  and  that  used  by  Lade  (1972)  in 
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triaxial  and  plane  strain  tests  are  shown  in  Figure  7.5. 
Sands  used  in  both  occasions  are  not  only  from  the  same 
source,  but  as  it  can  be  seen  in  that  Figure,  almost  have 
the  same  grain  size  distribution.  A  summary  of  the 
characteristics  of  these  sands  as  provided  by  Wong  (1978)  is 
given  in  Table  7.1. 

Separate  from  the  others  two  more  passive  pressure 
tests  were  carried  out  by  Wong  using  sand  which  contained 
thin  layers  of  dyed-red  sand.  These  tests  were  used  to 
identify  the  location  and  shapes  of  slip  planes.  Figure  7.6 
illustrates  the  observed  slip  planes  for  Dr=90%  at  the  end 
of  a  test  after  a  10  degree  wall  rotation.  What  is  shown  in 
this  Figure  is  the  sum  of  all  slip  planes  developed 
progress i ve 1 y  from  top  to  bottom  for  increasing  wall 
rotation.  For  example,  if  the  tests  were  carried  out  only  up 
to  5  degrees  rather  than  10  degrees  of  rotation,  the  lower 
slip  planes  would  not  yet  have  developed. 

Figure  7.7  and  Figure  7.8  show  average  normal  and  shear 
stresses  on  the  wall  as  well  as  the  coefficient  of  passive 
earth  pressure,  Kp,  for  two  different  densities  of  sand.  The 
distribution  of  these  stresses  along  the  wall  are  given  in 
Figure  7.9  and  Figure  7.10.  The  mobilization  of  stresses 
follow  a  pattern  which  can  be  summarized  as  follows. 
Stresses  acting  on  the  top  panel  of  the  wall  reach  to  a 
maximum  first.  With  increasing  wall  rotation,  failure 
spreads  along  the  wall  downward,  and  larger  portions  of  the 
soil  mass  away  from  the  wall  are  affected.  The  movement 
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below  the  mid  height  of  the  wall  is  not  sufficient  to 
develop  a  failure  state  in  the  soil.  As  shown  in  Figure  7.9 
and  Figure  7.10,  normal  stresses  on  the  third  and  fourth 
panels  do  not  show  any  sign  of  reaching  a  peak  even  at  8 
degrees  of  wall  rotation. 

The  requirements  for  strain  softening  to  take  place  are 
discussed  in  Chapter  5.  In  laboratory  tests,  for  example  see 
Figure  7.15,  dense  sands  under  low  confining  pressures  show 
significant  softening  response.  These  conditions  are 
satisfied  for  Wong's  earth  pressure  tests  when  dense  sand  is 
used.  Opposite  to  anticipations  though,  the  strain  softening 
does  not  seem  to  play  a  significant  role  on  the  measured 
normal  stresses  acting  on  the  wall  panels  as  shown  in  Figure 
7.10.  One  would  expect  a  larger  drop  in  stresses  beyond  the 
peak  stress  level.  Any  explanation  for  this  interesting 
observation  will  not  be  attempted  here,  but  it  is  left  to 
researchers  who  would  like  to  work  on  the  subject. 

In  any  analysis  based  on  rigid-plastic  soil  model,  it 
is  implied  that  an  infinitesimal  movement  along  a  failure 
surface  is  sufficient  to  develop  the  full  shear  strength. 
Consequently,  failure  is  reached  simultaneously  at  every 
point  within  the  deforming  part  of  a  soil  mass.  In  most 
analysis,  the  same  assumptions  are  also  made  for  the 
soil-wall  interface  behavior.  No  attention  is  paid  to  the 
magnitude  of  deformations.  A  suitable  factor  of  safety  is 
expected  to  prevent  excessive  movements.  For  the  analysis  of 
a  retaining  structure  subject  to  passive  earth  pressure,  it 
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is  conventional  to  draw  a  rupture  surface  that  passes 
through  the  toe  of  the  wall  and  extends  to  the  free  ground 
surface.  In  textbooks,  most  of  the  discussion  goes  around 
the  shape  of  this  assumed  rupture  surface.  The  observations 
by  Wong  (1978)  as  well  as  the  experiments  by  James  and 
Bransby  (1970)  contradict  the  basic  assumptions  of  limit 
equilibrium  method  of  analysis  for  this  mode  of  wall 
movement.  In  reality,  rupture  surfaces  do  not  start  from  the 
toe  of  a  wall  rotating  about  its  base.  Even  at  large 
deformations,  the  failure  zone  extends  only  from  the  surface 
to  somewhere  around  the  middle  of  the  wall.  This  is  a  good 
example  of  a  case  where  part  of  the  soil  mass  is  in  a  state 
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requires  both  the  deformation  properties  and  failure 
conditions  of  soils  to  be  represented  reasonably  well.  Based 
on  these  circumstances,  the  choice  of  the  passive  earth 
pressure  problem  for  the  evaluation  of  Lade's  model  seems  to 
be  a  proper  one. 


7.4  ANALYTICAL  WORK  BY  WONG 

In  the  analytical  part  of  his  work,  Wong  (1978)  made 
predictions  for  the  behavior  of  a  model  wall  by  using  four 
constitutive  models  which  were  considered  to  have  potential 
practical  value.  They  are: 

1)The  nonlinear  elastic  hyperbolic  model 
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2)  The  Mohr -Cou lomb  simple  el asto-pl ast ic  model 

3) The  s tress-di 1 atancy  model 

4)  Lade's  work  hardening  model 

A  summary  of  Wong's  predictions  in  relation  to  all  four 
models  is  reproduced  in  Figure  7.11  and  Figure  7.12  in  order 
to  compare  with  the  predictions  of  present  study.  As 
demonstrated  in  these  Figures,  the  wall  behavior  is 
predicted  quite  well  by  the  simple  elasto-plastic  model. 
Therefore,  a  brief  presentation  of  the  character i st i cs  of 
simple  elasto-plastic  model,  as  used  in  Wong's  analysis,  is 
considered  necessary.  The  model  employs  the  Mohr -Cou lomb 
yield  criterion  without  any  work  hardening.  Deformations  are 
elastic  up  to  failure  and  elasto-plastic  thereafter.  The 
flow  rule  used  is  non-associ ated .  Hence,  the  cost  of  running 
a  finite  element  analysis  is  very  similar  to  that  of  Lade's 
model.  The  secant  modulus  at  70%  of  the  maximum  deviatoric 
stress  is  used  as  the  elastic  modulus.  The  friction  angle  at 
peak  strength  as  well  as  the  elastic  modulus  are  obtained 
from  conventional  triaxial  compression  test  data.  It  should 
be  emphasized  here  that  the  finite  element  analysis  of  the 
experimental  wall  is  a  plane  strain  problem  and  the 
stress-strain  behavior  of  soils  in  conventional  triaxial 
tests  is  significantly  different  than  the  soil  behavior 
under  plane  strain  test  conditions. 

The  details  of  the  nonlinear  elastic  hyperbolic  model 
and  the  stress-di 1 atancy  model  will  not  be  given  here.  These 
and  the  soil  parameters  used  for  each  model  can  be  found  in 
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Wong  (  1978) . 

7.5  PREDICTIONS  BY  THE  FINITE  ELEMENT  ANALYSIS  OF  PRESENT 
STUDY 

A-PRELIMINARIES 

Before  a  proper  analysis  of  model  wall  behavior  was 
attempted  with  a  sufficient  number  of  elements,  the  finite 
element  program  of  the  present  study  was  used  to  predict  the 
coefficient  of  earth  pressure  versus  horizontal  strain 
relationship  similar  to  the  one  given  in  Figure  7.1.  This  is 
done  by  using  a  single  element.  The  selected  mode  of  wall 
movement  considered  here  is  a  pure  translation  rather  than  a 
rotation.  A  vertical  wall  with  frictionless  surface  is 
advanced  horizontally  to  induce  Rankine  states  of  failure  in 
soil.  Using  Lade's  work  hardening  model,  the  predictions  of 
the  strains  developed  for  the  stress  conditions  ranging  from 
Kc  at  rest  to  failure,  in  both  active  and  passive  sense, 
are  calculated.  The  results  are  shown  in  Figure  7.13.  The 
relative  orders  of  magnitude  of  earth  pressure  and  strain 
values  are  in  accordance  with  experimental  results  given  in 
the  literature.  Roscoe  and  Bur  land  (1968)  presented  similar 
predictions  for  wet  clays  to  illustrate  and  emphasize  the 
practical  significance  of  the  Cambridge  model.  However, 
Morgenstern  and  Eisenstein  (1970)  have  indicated  that  these 
relations  have  not  been  used  in  practice  as  yet.  It  should 
be  noted  that  the  relation  between  Kp  and  horizontal  strain 


■ 


1 1 

V 

s 

' 


* 


128 


is  a  function  of  the  state  of  stress  before  the  wall 
movement  is  initiated.  In  other  words,  the  passive  pressure 
coefficient  Kp  versus  horizontal  strain  €horz  relation 
varies  with  depth.  As  shown  in  Figure  6.7,  horizontal 
strains  to  failure  increase  with  increasing  depth.  It  is 
well  Known  that  the  strains  required  for  soil  samples  to 
reach  failure  increase  with  increasing  confining  pressures 
in  triaxial  tests  as  shown  in  Figure  4.8  and  Figure  4.9.  It 
appears  that,  this  character i st i cs  of  soil  behavior  is  also 
applicable  for  stress  paths  following  stress  states  in 
passive  earth  pressure  conditions.  However,  the  curves  given 
in  Figure  7.13  and  the  Kp versus  €h0rzrel at ionship  offered  by 
Roscoe  and  Burland  (1968)  are  calculated  for  highly 
idealized  boundary  conditions.  Even  under  these  hypothetical 
conditions  they  only  correspond  to  soil  behavior  at  one 
point  in  the  whole  soil  mass.  Therefore,  they  should  not  be 
expected  to  be  useful  in  great  detail  to  practicing 
engineers.  Lambe  and  Whitman  (1969)  use  triaxial  test  data 
for  the  same  problem  to  estimate  the  magnitude  of  strains. 
In  order  to  find  out  possible  consequences  of  using  triaxial 
data  for  the  analysis  of  passive  earth  pressure  problems, 
the  following  calculations  are  carried  out  by  using  Lade's 
work  hardening  model.  It  is  assumed  that  three  soil  samples 
are  obtained  from  locations  shown  in  Figure  6.6.  It  is  also 
assumed  that  K0  conditions  prevail  in  these  soil  elements, 
and  this  provides  the  initial  state  of  stress  before  the 
loading  commences.  During  subsequent  loading  vertical  stress 
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is  Kept  constant  but  all  around  horizontal  normal  stresses 
are  increased  until  failure  occured.  These  calculations  are 
carried  out  for  the  same  soil  used  for  the  passive  pressure 
calculations  which  are  shown  in  Figure  6.7.  The  results  of 
such  calculations  are  given  in  Figure  7.14.  A  comparison  of 
Figure  6.7  with  Figure  7.14  indicates  that  horizontal 
strains  required  to  reach  failure  are  about  50-100%  more  for 
plane  strain  calculations  than  triaxial  results.  Before  the 
failure  state  is  reached,  the  stress-strai n  curves  of 
triaxial  tests  have  steeper  slopes.  The  implication  of  this 
is  that  the  use  of  triaxial  data,  ( 0~ver  t  =const.  and  (7horz 
increasing  to  failure),  in  the  analysis  of  passive  earth 
pressure  problems  gives  deformation  results  which  are  on  the 
unsafe  side. 

B-ANALYSES  OF  THE  WALL  IN  PRESENT  STUDY 

In  the  present  study  two  analyses  were  performed  to 
predict  the  experimental  results  of  passive  earth  pressure 
tests  by  Wong  (1978).  The  predictions  were  made  for  both 
dense  and  medium  dense  sands. 

I  -  Ana lysi s  of  Wa 1 1  Behavior  for  Dense  Sand 

a-Soi 1  Parameters 

The  representative  soil  parameters  for  the  analysis  of 
wall  behavior  are  obtained  from  the  triaxial  test  data  which 
are  provided  by  Wong  (1978).  The  procedures  followed  to 
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calculate  soil  parameters  are  the  same  as  described  by  Lade 
and  Duncan  (1975).  Their  values  are  listed  in  Table  7.2. 
Because  Wong  employed  a  different  method  of  evaluation  than 
Lade  and  Duncan ( 1 975 ),( compare  the  parameters  in  Table  7.1 
and  Table  7.2),  the  parameters  of  the  present  study  are  not 
exactly  the  same  as  the  ones  used  by  Wong  (1978).  Using  the 
data  given  in  Table  7.2,  plane  strain  behavior  of  a  soil 
element  is  predicted  and  compared  with  experimental  results 
to  demonstrate  that  the  parameters  used  in  the  analyses  of 
present  study  are  appropriate.  This  comparison  is  shown  in 
F i gure  7.15. 
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b-Boundary  Conditions 

The  finite  element  mesh  and  the  boundary  conditions  for 
the  analysis  of  the  rigid  wall  are  shown  in  Figure  7.16.  For 
each  load  increment,  the  horizontal  displacements  of  the 
nodes  along  the  rotating  wall  and  soil  interface  are 
specified.  The  measurements  of  shear  stresses  given  in 
Figure  7.9  are  used  to  calculate  the  equivalent  shear  forces 
to  be  applied  to  the  nodes.  Thus,  the  wall  is  no  longer 
needed  in  the  analysis  of  this  problem. 

c-  Results  of  Analysis 

Results  of  the  analysis  involving  dense  sand  are  shown 
in  Figure  7.17  and  Figure  7.18.  In  these  figures  the  average 
normal  stress  as  well  as  its  distribution  along  the  wall  are 
plotted  against  the  angle  of  wall  rotation.  The  experimental 
results  and  the  predictions  made  by  Wong  are  also  plotted  on 
the  same  figures  for  the  purpose  of  comparison.  Additional 
curves  for  the  average  values  of  coefficients  of  passive 
earth  pressure  are  given  in  Figure  7.17..  If  the  variation  of 
Kp  on  each  panel  is  calculated  separately,  the  following 
numbers  would  be  obtained  according  to  the  present  study. 

Kp  va 1 ues  on  each  panel : 

(Note:The  numbers  in  parentheses  indicate  the  experimental 
values . ) 

At  0.5  degree  wall  rotation: 


PANEL  1  =  13  (13) 
PANEL  3=5  (4.6) 


PANEL  2=12(11) 
PANEL  4  =  0.9  (1.0) 
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At  1.0  degree  wall  rotation: 

PANEL  1  =  15  (15)  PANEL  2=15(14) 

PANEL  3  =  9  (7.4)  PANEL  4  =  1.9  (1.2) 

As  the  angle  of  wall  rotation  increases,  Kp  values  increase 
as  well,  but  not  at  the  same  rate  for  all  panels. 

The  growth  of  highly  stressed  zones  is  shown  in  Figure 

7.19. 

1 1  -  Ana  1 ys i s  of  Wa 1 1  Behavior  for  Medium  Dense  Sand 

a-  Soi 1  Parameters 

The  soil  parameters  used  for  medium  dense  sand  in  the 
analysis  of  the  wall  are  different  than  the  ones  provided  by 
Wong  (1978).  The  reasons  for  not  using  Wong's  parameters  are 
discussed  next. 

1- Wong  used  a  different  method  of  evaluation  for  the 
soil  parameters  than  Lade  and  Duncan  (1975). 

2- Wong' s  description  of  the  state  of  medium  dense  sand 
in  the  soil  bin  indicates  that  the  density  of  each  layer 
varied  within  its  depth.  Dense  and  loose  sublayers  were 
present  within  each  layer. 

3- There  was  another  difficulty  in  obtaining  a 
homogeneous  soil  deposit  in  medium  dense  sand.  This  is 
related  to  the  striations  as  photographed  by  Wong. 

4- Measurements  on  the  experimental  wall  by  Wong  showed 
that  the  coefficient  of  earth  pressure  at  rest  K0  was 
0.35  at  zero  wall  rotation.  This  value  of  coefficient  of 
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earth  pressure  corresponds  to  a  maximum  stress  level 
K-j  =64  rather  than  68.4  which  is  derived  from 
interpolated  triaxial  curves  by  Wong  (1978). 

Due  to  the  uncertainties  about  the  state  of  soil  in 
medium  dense  sand,  the  following  approach  is  taken  to 
evaluate  soil  parameters.  First  K-j  is  chosen  64.  Then  the 
stress-strain  curves  corresponding  to  K-)  are  obtained  by 
interpolating  between  the  triaxial  test  data  for  dense  and 
loose  sands.  After  that  the  procedures  given  in  Appendix  A 
are  followed  to  determine  the  soil  parameters.  The  results 
are  listed  in  Table  7.2. 

b-  Boundary  Conditions 

The  finite  element  mesh  for  the  analysis  of  the  wall  in 
medium  dense  sand  is  the  same  as  for  dense  sand  as  shown  in 
Figure  7.16.  Boundary  conditions  along  the  wall  are 
simulated  by  specifying  horizontal  displacements  and  the 
shear  forces  which  are  obtained  from  experimental  results 
given  i n  F igure  7.10. 

c-  Results  of  the  Analysis 

Figure  7.20  and  Figure  7.21  show  the  average  normal 
stress  and  its  distribution  along  the  wall  respectively.  The 
experimental  results  and  the  predictions  of  Wong  are  given 
on  the  same  figures.  The  average  value  of  the  coefficient  of 
passive  earth  pressure  is  illustrated  in  Figure  7.20.  The 
breakdown  of  Kp  into  each  panel  is  given  next. 
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Kp  va 1 ues  on  each  pane  1 : 

(Note:The  numbers  in  parentheses  indicate  the  experimental 
va 1 ues . ) 


At  0.5  degree  wa 1 1 
PANEL  1  =  7.3  (7.3) 

PANEL  3=3(3) 

At  1.0  degree  wa 1 1 
PANEL  1  =  8.5  (8.7) 

PANEL  3  =  4.6  (4.4) 


rotat i on : 

PANEL  2  =  5.8  (4.2) 

PANEL  4  =  0.7  (0.45) 

rotat i on : 

PANEL  2  =  7.5  (5.9) 

PANEL  4  =  1.4  (0.6) 


Growth  of  highly  stressed  zones  for  medium  dense  sand 
is  shown  in  Figure  7.22.  Contours  of  cumulative  shear  strain 
and  volumetric  strains  are  given  in  Figure  7.23. 


7.6  INTERPRETATION  OF  FINITE  ELEMENT  PREDICTIONS 

In  this  Section  particular  details  of  the  finite 
element  predictions  are  examined  in  the  following  order. 

I) The  significance  of  the  results  presented  in 
previous  sections  are  discussed  on  the  bases 
of  comparisions  between: 

a)  the  experimental  results  and  the  results  of 
present  study, 

b) the  calculations  of  Wong  and  the  predictions 
of  present  study. 

II) Several  important  factors  which  have  possible 

effects  on  the  results  are  elaborated. 

III ) Final ly,  in  the  light  of  evidence  presented 
above,  the  theoretical  and  practical  value  of 
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Lade's  model  is  discussed. 


I ) Comparative 

studies  of  average  normal  stresses 

on 

the 

i ns t rumen ted 

wall,  as  shown  i n 

F igure  7.17 

and  Figure  7 

.20, 

indicate  that 

the  predictions  of 

the  present 

study 

at 

al  1 

stages  of 

tests  are  within 

reasonable 

limits 

of 

the 

experimental 

data  for  both  dense 

and  medium 

dense  sands. 

The 

analysis  of  the  present  study  gave  average  stress  values 
which  were  not  different  from  the  experimental  results  by 
more  than  15-20%  for  any  specified  angle  of  wall  rotation. 
For  similar  conditions  Wong's  calculations  gave  average 
normal  stresses  reaching  up  to  1.7  to  1.9  times  the 
experimental  values.  Beyond  approximately  2.3  degrees  of 
wall  rotation  Wong  was  not  able  to  continue  his  finite 
element  analysis  for  numerical  reasons.  Calculations  of  the 
present  study  were  not  continued  any  further  than  2.2 
degrees  and  3.7  degrees  for  dense  and  medium  dense  sands 
respectively,  because  the  trend  had  already  been 
establ i shed . 

On  the  other  hand,  if  the  angle  of  wall  rotation  is  to 
be  predicted  for  a  specified  average  normal  stress,  the 
following  examples  will  be  more  relevant.  For  dense  sand,  at 
a  normal  stress  of  500  psf  ( =0 . 245kg/sq . cm) ,  the 
experimental  results  are  3.16  times  larger  than  the 
calculations  of  Wong  and  1.67  times  larger  than  the 
predictions  of  the  present  study.  For  medium  dense  sand  at  a 
normal  stress  of  300  psf  ( =0 . 1 48Kg/sq . cm ) ,  the  experimental 
results  are  4.01  times  larger  than  Wong's  calculations  and 
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1.6  times  larger  than  the  predictions  of  present  study. 
These  numbers  can  be  taken  as  an  indication  of: 

1) The  finite  element  analysis  of  the  present  study 
gave  much  better  predictions  than  Wong's 
ca leu  1  at  ions . 

2) The  predictions  of  the  present  study  are 
reasonably  good  for  both  low  and  high  stress  levels. 

3 )  A 1  though  the  results  of  the  present  study  are  very 
reasonable,  they  are  still  on  the  unsafe  side  of  the 
experimental  results. 

When  attention  is  turned  to  the  distribution  of  normal 
stresses  along  the  wall  a  better  understanding  of  the 
details  of  the  problem,  such  as  the  development  of  contact 
pressures,  can  be  obtained.  Experimental  results  indicate 
that  the  soil  elements  at  shallow  depths  behind  the  wall 
reach  the  maximum  stress  level  first.  As  the  angle  of 
rotation  increases,  the  failure  grows  deeper.  For  example, 
as  shown  in  Figure  7.18  and  Figure  7.21  the  soil  elements 
next  to  the  top  two  panels  reach  to  the  maximum  stress  level 
first  while  the  soil  against  lower  two  panels  is  still 
loading.  This  phenomenon  has  been  modelled  well  by  Lade's 
stress-s trai n  law  as  the  results  of  the  present  study  show. 
Yet  it  is  not  possible  to  observe  this  feature  in  Wong's 
calculations.  For  example,  in  his  calculations  the  soil  next 
to  Panel  1  in  both  dense  and  medium  dense  sands  does  not 
seem  to  be  getting  closer  to  failure  state.  Differences 
between  the  predictions  of  the  present  study  and  that  of 


’ 

' 


* 


. 

:  '  W  ~ 

' 


137 


Wong  are  most  likely  due  to  different  procedures  in  using 
the  hardening  law  of  the  model  in  these  analyses  (see 
Chapter  3 ) . 

Figure  7.19  and  Figure  7.22  show  the  expansion  of 
highly  stressed  zones  in  soil  mass  as  predicted  by  the 
present  study.  It  appears  that  highly  stressed  zones  in 
dense  sand  spread  out  more  than  that  of  medium  dense  sand 
for  the  same  degree  of  wall  rotation.  For  example  at  2 
degrees  of  rotation,  the  region  of  dense  sand  stressed  above 
fp  >  0.9 8  is  about  2  to  3  times  more  than  the 
cor respondi ng  area  for  medium  dense  sand.  It  should  be  noted 
that  the  gradual  growth  of  the  zone  of  failure  is  well 
represented  in  this  analysis.  Failure  has  not  been  reached 
in  all  soil  elements  (see  for  example  Figure  7.19)  along  the 
wall  face,  but  only  the  sand  within  one  third  of  the  wall 
height  from  the  surface  has  experienced  failure. 

Figure  7.23  shows  predictions  of  present  study  for  the 
cumulative  maximum  shear  and  volumetric  strains  at  the  end 
of  2  degrees  of  wall  rotation.  When  compared  with  curves 
measured  by  James  and  Bransby  (1970), (see  Figure  7.2),  the 
shear  strain  contours  given  in  Figure  7.23  indicate  that  the 
strains  within  the  soil  mass  are  predicted  very  well  by 
Lade's  model  for  complex  loading  conditions  . 

The  finite  element  analysis  of  the  wall  by  Wong  employs 
several  stress-strain  laws.  The  level  of  success  of  each 
model  is  shown  in  Figure  7.11  and  Figure  7.12.  If  these 
results  are  compared  with  the  predictions  of  the  present 
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study  which  are  given  in  Figure  7.18  and  Figure  7.21,  one 
can  conclude  that,  Lade's  model  gives  much  better  results 
than  the  hyperbolic  stress-strai n  law  and  the 
stress-di 1 atancy  model.  Comparisons  between  the  simple 
el asto-pl ast i c  model  and  Lade's  model  indicate  that  the 
predictions  related  to  the  top  three  panels  are  in  favor  of 
Lade' s  model . 

II)Factors  which  have  possible  effects  on  the  finite 
element  results  are  discussed  next.  These  are  related  to: 

a) Accuracy  of  experimental  results 

b)  Influence  of  soil  parameters 

c )  Approximat ions  related  to  the  finite  element 

ana  lysis. 

a ) Accuracy  of  exper i menta 1  data 

As  will  be  quoted  next  from  Wong  (1978),  there  were 
several  difficulties  related  to  the  experimental 
work . 

"Good  uniformity  of  sand  density  was  obtained  except  in 
the  areas  adjacent  to  the  wall... 

Although  the  measured  densities  appeared  to  be  uniform 
throughout  the  bin,  the  density  of  each  thin  layer  varied 
within  its  depth.  The  loose  and  dense  sublayers  were  more 
pronounced  in  loose  samples  than  in  dense  samples... 

Each  layer  deposited  in  the  medium  dense  samples 
(Dr=65%)  was  about  0.3in.  thick,  and  those  in  the  dense 
( D r  =  90% )  samples  were  about  0.07in.  thick.  Striations  were 
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readily  identifiable  within  the  medium  dense  sample,..." 

In  addition  to  the  difficulties  involved  in  achiving  a 
uniform  sand  deposit,  there  was  another  problem  related  to 
the  rigidity  of  the  testing  apparatus.  This  will  be  quoted 
again  from  Wong  (1978). 

"It  was  found  that  the  bottom  of  the  wall  tended  to 
move  away  from  the  sand  during  rotation;  the  amount, 
however,  was  very  small.  When  the  tank  was  filled  with  dense 
Monterey  No.O  sand,  the  movement  was  less  than  0.03in..." 

Although  it  is  quite  possible  to  model  such  a  wall 
movement  in  the  finite  element  analysis,  this  was  not 
attempted  in  the  present  study.  Yet  its  importance  may  not 
be  ignored.  If  this  movement  was  not  allowed  in  experiments, 
the  bottom  panel  of  the  wall  would  have  carried  larger  load 
than  what  is  shown  in  Figure  7.18  and  Figure  7.21.  This 
would  allow  the  predictions  of  the  present  study  to  be 
closer  to  the  experimental  results. 

b ) Inf luence  of  soi 1  parameters 

The  parameters  used  for  soil  have  major  influence  on 
the  outcome  of  the  finite  element  analysis.  It  should  be 
noted  that  the  selection  of  the  parameters  for  Lade's  model 
is  not  based  on  the  past  experience;  they  have  to  be 
obtained  from  conventional  triaxial  test  data  by  following 
the  steps  given  in  Appendix  A.  In  soil  testing  repeatability 
of  tests  within  10%  error  is  considered  to  be  normal.  If  the 
triaxial  compression  test  data  given  in  Figure  A.1  in 
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Appendix  A  is  examined,  one  can  observe  that  measured 
maximum  stress  level  Ki  for  three  tests  with  different 
confining  pressures  vary  from  K1 =78.67  for  (73=0.6  Kg/sq.cm 
to  K-j  =88 . 9  for  (73=0.3  Kg/sq.cm.  For  confining  pressure  1.2 
Kg/sq.cm,  K-j  is  80.95.  The  parameter  K-j  used  in  the  present 
analysis  is  the  average  of  all  three  K1  values  which  is 
83.6.  If  Ki  obtained  from  triaxial  test  with 

0*3  =0 . 3kg/sq .  cm .  was  excluded,  a  lower  value  would  be  used 
for  K1  and  that  would  bring  predictions  of  present  study 
closer  to  experimental  results. 

c ) Approximations  related  to  the  finite  element  ana  lysis 

It  is  most  likely  that  a  finer  mesh  would  further 
increase  the  flexibility  of  soil  mass  in  the  finite  element 
analysis.  Especially,  where  the  gradient  of  stress  change  is 
large  ,  an  increased  number  of  elements  would  improve  the 
predictions.  Such  a  point  in  the  analysis  of  the  wall  is 
around  the  toe. 

Finally,  the  procedure  used  for  the  nonlinear  analysis 
may  have  an  effect  on  the  results.  It  is  well  Known  that  the 
tangent  stiffness  method  gives  results  which  are  stiffer 
than  the  correct  solution. 

It  is  believed  that  these  factors  listed  above  have 
some  influence  on  the  predictions  of  present  study.  The 
investigation  of  their  relative  importance  requires  more 
effort.  Subject  to  the  problems  mentioned  above,  conclusions 
will  be  drawn  next  for  the  value  of  Lade's  work  hardening 
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mode  1 . 

1 1 1 ) Eva  1 uat i on  of  the  model 

The  main  objective  of  the  work  presented  in  this 
Chapter  is  to  investigate  the  degree  of  agreement  between 
the  predictions  made  by  Lade's  model  and  the  experimental 
results  for  complex  problems.  Comparisons  indicate  a  very 
good  agreement  in  general  between  the  observed  stress 
distribution  on  the  wall  and  the  measured  values. 
Predictions  of  zones  of  high  stress  level  and  their  growth 
as  a  function  of  wall  rotation  agree  well  with  the  measured 
values.  Predictions  at  working  load  range  as  well  as  at  high 
stress  levels  are  rather  impressive.  It  seems  that  the 
success  in  predictions  is  due  to  the  abilities  of  the  model 
to  predict  the  path  dependent  soil  behavior  and  the  failure 
conditions  reasonably  well.  These  results  show  that  a 
unified  analysis  of  deformations  and  stability  can 
successfully  be  carried  out.  To  conduct  a  reliable  analysis 
though,  the  soil  model  used  has  to  be  a  capable  one. 
Comparisons  of  predictions  made  by  the  model  used  here  and 
the  other  models  selected  by  Wong  show  clearly  that  Lade's 
model  gives  better  results  than  simple  elasto-plastic, 
nonlinear  elastic  and  st ress-di 1 atancy  models.  In  this 
conclusion  it  is  assumed  that  Wong's  calculations  for  these 


models  are  correct. 
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Figure  7.1.  Relation  between  wall  rotation  and  earth 
pressure  for  different  densities  (After  Terzaghi ,  1954). 
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0  :  Angle  of  wall  rotation  about  the  toe 


Figure  7.3.  Rupture  planes  observed  on  radiographs  exposed 
at  three  stages  of  wall  rotation  (After  James  and  Bransby, 
1970) . 
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Figure  7.5.  Gradation  curves  for  Monterey  No.O  sand  (After 
Wong ,  1978). 


TABLE  7. 1 


STRESS-STRAIN  PARAMETERS  OBTAINED  BY  WONG  (1978) 

FOR  MONTEREY  No.O  SAND. 


Parameter 

Dense 

Sand 

Med.  Dense 
Sand 

Kur 

3720 

2330 

n 

0.78 

0.79 

V 

0.2 

0.26 

K, 

83.6 

68.4 

ft 

30.3 

30.9 

A, 

0.57  2 

0.451 

a2 

7.11 

12.96 

l 

1 .45 

1.38 

M 

0.0000746 

0.000196 

rf 

0.97 

0.99 
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Figure  7.6.  Observed  rupture  planes  at  10  degree  wall 
rotation  (After  Wong,  1978). 
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Figure  7.7  Experimental  results  for  dense  sand.  Variations 
of  average  normal  and  shear  stresses  and  passive  earth 
pressure  coefficient  with  the  angle  of  wall  rotation.  (After 
Wong ,  1978). 
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Figure  7.8  Experimental  results  for  medium  dense  sand. 
Variations  of  average  normal  and  shear  stresses  and  passive 
earth  pressure  coefficient  with  the  angle  of  wall  rotation. 

( After  Wong ,  1978). 
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Figure  7.9  Experimental  results  for  dense  sand.  Distribution 
of  normal  and  shear  stresses  on  four  panels  of  the  wall. 

( After  Wong ,  1978 ) . 
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7.10  Experimental  results  for  medium  dense  sand, 
bution  of  normal  and  shear  stresses  on  four  panels  of 


the  wall.  (After  Wong,  1978). 
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Figure  7.13  Relation  between  active  and  passive  earth 
pressure  coefficient  and  horizontal  strain. 
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Figure  7.14  Predictions  for  (Tvert  =const. tests  with  C7horz. 


increasing  to  failure. 
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TABLE  7.2 


STRESS-STRAIN  PARAMETERS  USED 


IN  PRESENT  STUDY 


FOR  MONTEREY  No.O  SAND 


Parameter 

Dense 

Sand 

Med.  Dense 
Sand 

Kur 

2800 

169  3 

n 

0.  8 

0.85 

V 

0 

0 

Ki 

83.6 

64.0 

ft 

38 

33.9 

Ai 

0.45 

0.397 

A2 

14.85 

16.28 

l 

1.686 

1  .19 

M 

0.00025 

0.00062 

rt 

0.973 

0.968 

(%)  ( (7,  -  CT3 )  kg/cm 
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Figure  7.15  Experimental  results  and  predictions  of  present 
study  for  dense  Monterey  No.O  sand  in  plane  strain  tests. 
(Experimental  data  after  Lade,  1972) 
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Figure  7.16  Finite  element  mesh  and  boundary  conditions  for 
the  analysis  of  rigid  wall. 
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ANGLE  OF  ROTATION  (Degree) 


Figure  7.17  Variation  of  average  normal  stress  with  wall 
rotation  for  dense  Monterey  No.O  sand. 
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Figure  7.19  Zones  of  high  stress  level  as  a  function  of  wall 
rotation  for  dense  Monterey  No.  0  sand. 
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Figure  7.20  Variation  of  average  normal  stress  with  wall 
rotation  for  medium  dense  Monterey  No.O  sand. 
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Figure  7.22  Zones  of  high  stress  level  as  a  function  of  wall 
rotation  for  medium  dense  Monterey  No.  0  sand. 
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dense  Monterey  No.  0  sand. 
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CHAPTER  8 


SUMMARY  AND  CONCLUSIONS 


8. 1  SUMMARY 

The  present  study  is  an  attempt  to  reassess  the 
theoretical  and  practical  value  of  Lade's  work  hardening 
constitutive  law.  It  is,  therefore,  relevant  to  summarize 
the  char acter i st i cs ,  capabilities  and  limitations  of  the 
model . 

CHARACTERISTICS  OF  WORK  HARDENING  MODEL 

Lade's  first  model  which  is  evaluated  in  the  present 
study,  is  an  e 1 asto-pl ast i c  work  hardening  stress-strai n 
relationship.  It  is  developed  for  cohesionless  soils. 
Essential  components  of  the  model  are  the  yield  and  failure 
criteria,  flow  rule  and  hardening  law.  The  yield  surface  has 
the  shape  of  a  cone  which  expands  with  increasing  stress 
levels  and  eventually  reaches  the  failure  surface.  The 
effect  of  intermediate  principal  stress  in  the  development 
of  state  of  failure  is  accounted  for.  A  non-associ ated  flow 
rule  is  used  for  better  modelling.  Hence,  the  plastic 
potential  and  the  yield  surface  are  not  identical.  The  work 
hardening  law  is  based  on  an  experimentally  determined 
relationship  between  plastic  work  and  stress  level.  Details 
of  the  formulation  of  the  model  can  be  found  in  Chapter  2  or 
in  the  original  work  by  Lade  (1972). 
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CAPABILITIES  OF  WORK  HARDENING  MODEL 

The  stress-strain  law  is  designed  for  modelling  the 
most  essential  features  of  soil  behavior  such  as 
nonlinearity,  inelasticity,  path  dependancy  and  shear 
dilatancy.  The  effect  of  intermediate  principal  stress  on 
soil  behavior  is  included  in  the  formulation.  Coincidence  of 
strain  increment  and  stress  increment  axes  at  low  stress 
levels  with  transition  to  coincidence  of  strain  increment 
and  stress  axes  at  high  stress  levels  is  another  significant 
feature  of  the  model.  The  soil  behavior  under  reorientation 
of  principal  stress  axes  can  be  predicted  for  stress  paths 
with  increasing  stress  levels.  Comparisons  between 
predictions  of  the  model  and  the  observations  in  true 
triaxial,  plane  strain,  torsion  shear  and  conventional 
triaxial  apparatus  (used  in  both  compression  and  tension) 
have  demonstrated  that  the  stress-strain  law  gives  good 
results  for  all  stages  of  loading  up  to  and  including 
f a i 1 ure . 

Once  the  capabilities  of  the  model  were  understood  and 
the  level  of  success  in  predictions  was  established,  the 
attention  is  focused  on  the  limitations  of  the  model. 


LIMITATIONS  OF  WORK  HARDENING  MODEL 

As  presented  in  Section  2.3  the  limitations  of  the 
stress-strain  law  are: 

1.  Inability  to  model  strain  softening  (conditions  for 
strain  softening  to  take  place  are  given  in  Chapter  4). 
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2.  Inexistence  of  a  cap  yield  surface  results  in  inadequate 
modelling  for  certain  stress  paths  (see  Chapter  5). 

3.  Yield  and  failure  surfaces  are  conical  rather  than  the 
shape  of  a  bullet.  A  better  fit  to  experimental  results 
may  be  obtained  with  bullet  shape  surfaces  (see  Chapter 
4)  . 


4. 


5. 


6. 


Poisson's  ratio  is  equal  to 
the  predictions  for  elastic 
correspond  with  experiments 
If  the  stress  level 
reorientation  of  principal 
cannot  be  modelled  accurate 
Under  cyclic  loading  within 
level,  the  predictions  for 
Experiments  show  that  pi 
large  stress  reversals. 


zero.  Due  to  this  limitation 
strain  increments  do  not 
1  resu 1 t s  very  well, 
remains  the  same  during 
stress  axes,  soil  behavior 


ly. 

the  bounds  of  a  given  stress 
strains  are  purely  elastic, 
astic  strains  take  place  for 


REMARKS  ON  PREVIOUS  FINITE  ELEMENT  ANALYSES 


Since  the  first 

development  of  the 

model , 

several 

modi fi cat  ions 

have 

been  proposed  to 

the 

original 

formu 1  at i on . 

These  are 

discussed  in  Chapter  3. 

For  the 

analysis  of 

boundary 

value  problems,  a 

finite 

e 1 ement 

program  was 

developed 

1  by  Ozawa  (1974). 

A 1  1 

previous 

modifications  to  the  model,  before  the  present  study,  were 
implemented  into  Ozawa's  finite  element  program  by  Wong 
(1978).  For  the  reasons  discussed  in  Chapter  3,  the  results 
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of  the  finite  element  analysis  presented  by  Ozawa  (1974)  and 
Wong  (1978)  are  not  correct.  Therefore,  their  conclusions 
about  the  significance  of  Lade's  work  hardening  model  are 
not  relevant . 


8.2  CONCLUSIONS  OF  PRESENT  STUDY 

During  the  course  of  present  work  various  conclusions 
have  been  reached  in  connection  with: 

1.  The  formulation  of  work  hardening  model 

2.  The  hardening  law  of  strain  softening  model 

3.  A  possible  improvement  in  work  hardening  model 

4.  Finite  element  analysis 

5.  Passive  earth  pressure  problem. 

These  conclusions  are  presented  next. 


FORMULATION  OF  WORK  HARDENING  MODEL 

In  order  to  find  reasons  for  unsatisfactory  results  of 
the  finite  element  analyses  of  the  past,  the  original 
formulation  of  the  model  and  the  previous  modifications  are 
re-evaluated.  It  is  found  that  the  work  hardening  law  is  the 
source  of  difficulties  of  the  past.  It  is  emphasized  here 
that  the  relationship  between  stress  level  and  plastic  work 
depends  on  the  confining  pressure  for  conventional  triaxial 
compression  tests.  For  a  general  three  dimensional  problem, 
a  unique  relation  between  stress  level  and  plastic  work  does 
not  exist.  To  make  use  of  the  f  versus  Wp  relation  obtained 
from  triaxial  tests  in  the  analysis  of  a  general  3-D 
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problem,  a  somewhat  arbitrary  assumption  has  to  be  made. 
This  assumption  simply  suggests  that,  for  a  general  stress 
change  from  a  Known  stress  level  and  O3  value,  the  slope 
of  the  f  versus  Wp  curve  of  a  conventional  triaxial  test  is 
to  be  used  to  relate  the  increment  of  stress  level  to  the 
plastic  work  increment.  In  triaxial  data,  at  the  point  where 
the  slope  of  the  f  versus  Wp  curve  is  considered,  the  value 
of  (J 3  and  f  have  to  correspond  to  that  of  the  soil 
element  under  general  state  of  stress  (see  Chapter  3).  Using 
that  assumption,  measured  stress-strain  behavior  of 
cohesionless  soils  can  be  successfully  predicted.  At  the 
same  time,  the  problems  that  are  created  in  the  finite 
element  analysis  can  be  avoided. 

It  is  also  shown  here  that  the  modification  suggested 
by  Ozawa  (1974)  in  relation  to  parameter  a  adversely  affects 
the  capabilities  of  the  model  for  certain  stress  paths  (see 
Chapter  3 ) . 

The  use  of  modified  plastic  potential  function  which  is 
proposed  by  Ozawa  (1974)  produces  results  identical  to  those 
obtained  from  Lade's  original  formulation.  Therefore,  it  is 
unnecessary . 

HARDENING  LAW  OF  STRAIN  SOFTENING  MODEL 

Hardening-softening  law  of  Lade's  strain  softening 
model  is  based  on  a  relationship  between  stress  level  and 
total  plastic  work.  This  relationship  is  not  unique  but  is  a 
function  of  the  confining  pressure  for  conventional  triaxial 
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compression  tests.  In  order  to  use  such  a  relation  for 
general  stress  changes  and  also  to  avoid  problems  in  finite 
element  analysis,  an  interpretation  of  the  law  is  proposed 
( see  Chapter  4  )  . 

POSSIBLE  IMPROVEMENT  IN  WORK  HARDENING  MODEL 

One  of  the  limitations  of  the  work  hardening  model  is 
unsatisfactory  modelling  for  certain  stress  paths  as 
discussed  in  Chapter  5.  To  eliminate  this  limitation,  a  cap 
type  yield  surface  is  incorporated  into  Lade's  work 
hardening  model.  The  improvement  that  can  be  obtained  in 
predictions  is  demonstrated  by  comparing  experimental  and 
calculated  results. 

FINITE  ELEMENT  ANALYSIS 

Related  to  the  finite  element  analysis,  the  conclusions 
reached  in  different  sections  of  the  present  study  follow. 

1.  It  has  already  been  emphasized  that  the  proper  use  of 
the  hardening  law  of  Lade's  work  hardening  model  is 
essential  for  the  results  of  finite  element  analysis  to 
be  reliable  (see  Chapter  3). 

2.  Constitutive  laws  dealing  with  several  stress-strai n 
relationships  within  the  model  may  be  formulated  in  a 
more  suitable  way  for  finite  element  programming  (see 
Chapter  5 ) . 

3.  A  method  to  test  the  accuracy  and  the  correctness  of 
newly  developed  finite  element  programs  is  proposed  in 
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Chapter  6.  This  is  a  better  technique  than  what  has  been 
used  in  geotechnical  practice. 

PASSIVE  EARTH  PRESSURE  PROBLEM 

Preliminary  calculations  have  shown  that  the  passive 
earth  pressure  coefficient  versus  horizontal  strain  relation 
is  a  function  of  the  state  of  stress  at  rest  conditions.  In 
other  words,  it  depends  strongly  on  the  depth  of  soil 
element  under  consideration. 

Also,  the  analysis  of  passive  earth  pressure  problems 
should  not  be  carried  out  on  the  basis  of  triaxial  tests 
where  soil  samples  are  subjected  to  constant  vertical  stress 
and  increasing  horizontal  stresses  until  failure  occurs.  The 
use  of  such  triaxial  test  data  gives  deformation  results 
which  are  on  the  unsafe  side. 

The  finite  element  analysis  of  the  instrumented  wall 
has  indicated  that  Lade's  work  hardening  model  can  be  used 
successfully  for  deformation  analyses.  Comparative  studies 
of  experimental  and  analytical  results  of  passive  earth 
pressure  tests  show  that  the  model  is  capable  of  predicting 
the  most  important  features  of  the  wall  behavior  such  as  the 
distribution  of  contact  pressure  and  the  magnitude  of 
deformations  reasonably  well. 

When  compared  with  other  models  such  as  the  hyperbolic 
stress-strai n  law  and  the  stress-di 1 atancy  model,  Lade's 
model  gave  much  better  results  for  the  wall  behavior. 
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8.3  CLOSING  REMARKS 

In  the  present  study,  the  first  time  successful 
application  of  Lade's  work  hardening  model  to  a  complex 
problem  has  been  achieved.  Undoubtedly,  the  model  works  well 
at  all  stages  of  loading.  Therefore,  it  has  the  potential  to 
be  very  useful  to  engineering  practice.  It  is  hoped  that  the 
model  will  attract  more  attention  in  the  future. 


- 
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APPENDIX  A 


DERIVATION  OF  SOIL  PARAMETERS 

The  procedures  to  derive  stress-strain  parameters  of 
Lade's  work  hardening  model  are  presented  here  by  using 
dense  Monterey  No.O  sand  ( e  =  0 . 5  7  7 ) .  There  are  nine 
parameters  in  this  model,  and  they  can  be  obtained  from 
drained  conventional  triaxial  compression  test  data.  A 
minimum  of  two  and  preferably  three  or  more  triaxial  tests 
with  different  confining  pressures  are  required.  For  the 
evaluation  of  elastic  parameters,  a  minimum  of  two 
unloading-reloading  branches  are  necessary.  Figure  A-1  shows 
three  triaxial  test  data  which  are  used  to  derive  the 
parameters . 

1 .  DERIVATION  OF  PARAMETERS  FOR  ELASTIC  STRAINS 

The  procedures  described  here  for  elastic 
parameters  can  also  be  found  in  Duncan  and  Chang  (1970). 

Unloading-reloading  branches  of  the  curves  given  in 
Figure  A-1  are  used  first  to  calculate  the  elastic 
modulus,  Eur .  In  a  log-log  plot,  Eur  versus  (J3 
relation  is  shown  in  Figure  A-2.  The  slope  of  the 
straight  line  connecting  the  data  points  gives  the  value 
of  parameter  n.  The  modulus  value  at  confining  pressure 
equal  to  unity  provides  the  magnitude  of  Kur . 
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2.  DERIVATION  OF  PARAMETERS  RELATED  TO  PLASTIC  STRAINS 

There  are  six  parameters  for  the  calculation  of 
plastic  strain  increments.  Derivation  of  these 
parameters  follows. 

a ) PARAMETER  K1  IN  FAILURE  CRITERION 

Stress  level,  f,  corresponding  to  peak  strength  is 
calculated  for  each  stress-strain  curve  shown  in  Figure  A-1. 
The  results  are: 


K1 =81 .24 

for 

°3 

=  1  . 2kg/ sq . cm 

K1 =78.95 

for 

^3 

=0 . 6kg/sq . cm 

CO 

CD 

00 

00 

II 

for 

^3 

=0 . 3kg/ sq . cm 

The  mean  value 

of  K1 

i  s 

used  in  the  analysis  of 

instrumented  wall. 


b ) PARAMETER  A  IN  PLASTIC  POTENTIAL 


Parameter 

A  is  introduced  in 

Eqn 

2 . 6  for 

the 

calculation  of 

4-J 

l — 1 

CNI 

is  related  to 

the  direction 

of 

plastic  strain 

i ncrements . 

Steps  to 

be 

fol lowed  for 

the 

calculation  of 

A  are  given 

next . 

-Calculate  elastic  strain 

i ncrements 

( use 

Eqn .  2.3). 

-Subtract  elastic  strain  increments  from  total  strain 
increments  to  find  plastic  strain  increments. 

-For  each  stress  increment,  calculate  Vp  which  is  defined 
next . 


.  .  .  . Eqn.  A.  1 
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-The  next  step  requires  the  following  calculations. 
An  explicit  form  for  Vp  is  derived  first  by  using 
Eqn .2.5  and  Eqn .  2.7. 


-  V 


. . . . Eqn .  A . 2 


If  Eqn.  A. 2  is  rearranged,  an  expression  can  be  obtained 
for  K2. 


3.1, ( 1  +  v  ) 


K2  =  — ^ - — 

0*3  ( Oj  +  Oj-CJ3 ) 


Eqn .  A . 3 


All  variables  on  the  right  hand  side  of  Eqn.  A. 3  are  Known. 
-Plot  f<2  versus  f  for  each  load  increment.  As  shown  in 
Figure  A-3,  the  slope  of  K2  versus  f  relation  gives 
the  parameter  A. 

c ) PARAMETERS  RELATED  TO  WORK  HARDENING  LAW 

Experimental  results  by  Lade  (1972)  show  that  for  a 
constant  confining  pressure,  there  is  a  unique  relation 
between  Wp  and  f.  The  values  of  Wp  and  f  calculated  from 
experimental  data  for  dense  Monterey  No.O  sand  are  shown  in 
Figure  A. 4.  As  a  convenience  in  fitting  curves,  the  stress 
level  at  which  Wp  is  practically  equal  to  zero  is  called 
treshold  stress  level,  ft  .  As  shown  in  Figure  A.  4,  ft  is 
chosen  38.  The  relation  between  Wp  and  (f-ft  )  is 
approximated  by  hyperbolas,  for  which  the  following 


nt 


. 
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expression  is  used. 


(f-ft)=  - — —  •  •  •  Eqn  .  A.  4 

a  +  b.Wp 

The  values  of  a  and  b  can  be  derived  from  a  plot  on 
transformed  axes  as  shown  in  Figure.  A. 5.  In  this  Figure, 
is  the  intercept  of  a  straight  line  with  the  vertical  axis, 
and  b  is  the  slope  of  that  straight  line.  As  shown  in 
Figure.  2.3,  the  reciprocal  of  b  is  equal  to  ( f  -  f ^  ^  .  In 

order  to  find  the  parameter  rj.  ,  Eqn.  2.11  is  utilized;  An 
average  value  of  r^  is  used  in  the  finite  element  analysis 
of  rotating  wall  (see  Table  7.2). 

The  value  of  a  increases  with  increasing  confining 
pressures .  This  relation  is  given  by  Eqn.  2.9.  The 
parameters  M  and  1  of  Eqn.  2.9  are  derived  from  a  log- log 
plot  for  a  versus  (J3  relation.  Figure  A. 6  shows  the 
variation  of  a  with  CJ3  for  dense  Monterey  No.O  sand. 

Summary  of  all  nine  parameters  is  provided  in  Table 


7.2. 
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Figure  A.1  Three  drained  conventional  triaxial  compression 
tests  for  dense  Monterey  No.O  sand.  (After  Lade,  1972) 
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Figure  A. 2  Derivation  of  parameters  for  elastic  strains 
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Figure  A. 3  Variation  of 
Monterey  No.O  sand. 
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Figure  A. 4  Variation  of  plastic  work  for  dense  Monterey  No. 
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Figure  A. 5  Transformed  plastic  work  curves  for  dense 
Monterey  No.O  sand. 
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Figure  A. 6  Variation  of  initial  slope  of  plastic  curves  with 
confining  pressure  for  dense  Monterey  No.O  sand. 
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